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Chapter 1 

The Quantum Hall Effect 



1.1 Introduction 

The quantum Hall effect (QHE) is one of the most remarkable condensed-matter 
phenomena discovered in the second half of the 20th century. It rivals super- 
conductivity in its fundamental significance as a manifestation of quantum me- 
chanics on macroscopic scales. The basic experimental observation is the nearly 
vanishing dissipation 



of a real (as opposed to some theorist's fantasy) transistor-like device (similar in 
some cases to the transistors in computer chips) containing a two-dimensional 
electron gas subjected to a strong magnetic field. This quantization is universal 
and independent of all microscopic details such as the type of semiconductor 
material, the purity of the sample, the precise value of the magnetic field, and so 
forth. As a result, the effect is now used to maintain^ the standard of electrical 
resistance by metrology laboratories around the world. In addition, since the 
speed of light is now defined, a measurement of /h is equivalent to a mea- 
surement of the fine structure constant of fundamental importance in quantum 
elect. r()(l\nauii(:s. 

^Maintain docs not mean define. The SI ohm is defined in terms of the kilogram, the 
second and the speed of Ught (formerly the meter). It is best realized using the reactive 
impedance of a capacitor whose capacitance is computed from first principles. This is an 
extremely tedious procedure and the QHE is a very convenient method for realizing a fixed, 
reproducible impedance to check for drifts of resistance standards. It does not however define 
the ohm. Eq. (1.2) is given in cgs units. When converted to SI units the quantum of resistance 

is ^/e^(cgs) ~> ^ 25, 812.80 f2 (SI) where a is the fine structure constant and Z = y^//o7^ 
is the impedance of free space. 







(1.1) 



and the quantization of the Hall conductance 



a. 




(1.2) 
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In the so-called integer quantum Hall effect (IQHE) discovered by von Kl- 
itzing in 1980, the quantum number is a simple integer with a precision of 
about 10~^° and an absolute accuracy of about 10~^ (both being limited by our 
ability to do resistance metrology). 

In 1982, Tsui, Stormer and Gossard discovered that in certain devices with 
reduced (but still non-zero) disorder, the quantum number v could take on ra- 
tional fractional values. This so-called fractional quantum Hall effect (FQHE) 
is the result of quite different underlying physics involving strong Coulomb in- 
teractions and correlations among the electrons. The particles condense into 
special quantum states whose excitations have the bizarre property of being 
described by fractional quantum numbers, including fractional charge and frac- 
tional statistics that are intermediate between ordinary Bose and Fermi statis- 
tics. The FQHE has proven to be a rich and surprising arena for the testing 
of our understanding of strongly correlated quantum systems. With a simple 
twist of a dial on her apparatus, the quantum Hall experimentalist can cause 
the electrons to condense into a bewildering array of new 'vacua', each of which 
is described by a different quantum field theory. The novel order parameters 
describing each of these phases are completely unprecedented. 

Wc begin with a brief description of why two- dimensionality is important to 
the universality of the result and how modern semiconductor processing tech- 
niques can be used to generate a nearly ideal two-dimensional electron gas 
(2DEG). We then give a review of the classical and semi-classical theories of 
the motion of charged particles in a magnetic field. Next we consider the limit 
of low temperatures and strong fields where a full quantum treatment of the 
dynamics is required. After that wc will be in a position to understand the 
localization phase transition in the IQHE. We will then study the origins of 
the FQHE and the physics described by the novel wave function invented by 
Robert Laughlin to describe the special condensed state of the electrons. Finally 
we will discuss topological excitations and broken symmetries in quantum Hall 
ferromagnets. 

The review presented here is by no means complete. It is primarily an 
introduction to the basics followed by a more advanced discussion of recent 
developments in quantum Hall ferromagnetism. Among the many topics which 
receive little or no discussion are the FQHE hierarchical states, interlayer drag 
effects, FQHE edge state tunneling and the composite boson [1] and fermion [2] 
pictures of the FQHE. A number of general reviews exist which the reader may 
be interested in consulting [3-11] 

1.1.1 Why 2D Is Important 

As one learns in the study of scaling in the localization transition, resistivity 
(which is what theorists calculate) and resistance (which is what experimental- 
ists measure) for classical systems (in the shape of a hypercube) of size L are 
related by [12, 13] 



(1.3) 
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Figure 1.1: Schematic illustration of a GaAs/ AlAs heterostructure quantum 
well. The vertical axis is band energy and the horizontal axis is position in 
the MBE growth direction. The dark circles indicate the Si+ ions which have 
donated electrons into the quantum well. The lowest electric subband wave 
function of the quantum well is illustrated by the dashed line. It is common to 
use an alloy of GaAs and AlAs rather than pure AlAs for the barrier region as 
illustrated here. 



Two dimensions is therefore special since in this case the resistance of the sample 
is scale invariant and (e'^/h)R is dimcnsionlcss. This turns out to be crucial to 
the universality of the result. In particular it means that one does not have to 
measure the physical dimensions of the sample to one part in 10^° in order to 
obtain the resistivity to that precision. Since the locations of the edges of the 
sample are not well-defined enough to even contemplate such a measurement, 
this is a very fortunate feature of having available a 2DEG. It further turns 
out that, since the dissipation is nearly zero in the QHE states, even the shape 
of the sample and the precise location of the Hall voltage probes are almost 
completely irrelevant. 



1.1.2 Constructing the 2DEG 

There are a variety of technicjucs to construct two-dimensional electron gases. 
Fig. (1.1) shows one example in which the energy bands in a GaAs/AlAs het- 
erostructure are used to create a 'quantum well'. Electrons from a Si donor 
layer fall into the quantum well to create the 2DEG. The energy level ('electric 
subband') spacing for the 'particle in a box' states of the well can be of order 
10^ K which is much larger than the cryogenic temperatures at which QHE 
experiments are performed. Hence all the electrons arc frozen into the lowest 
electric subband (if this is consistent with the Pauli principle) but remain free 
to move in the plane of the GaAs layer forming the well. The dynamics of the 
electrons is therefore effectively two-dimensional even though the quantum well 
is not literally two-dimensional. 
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Heterostructures that are grown one atomic layer at a time by Molecular 
Beam Epitaxy (MBE) are nearly perfectly ordered on the atomic scale. In ad- 
dition the Si donor layer can be set back a considerable distance (~ 0.5/im) to 
minimize the random scattering from the ionized Si donors. Using these tech- 
niques, electron mobilities of 10^ cm^/Vs can be achieved at low temperatures 
corresponding to incredibly long mean free paths of ^ 0.1 mm. As a result of 
the extremely low disorder in these systems, subtle electronic correlation ener- 
gies come to the fore and yield a remarkable variety of quantum ground states, 
some of which we shall explore here. 

The same MBE and remote doping technology is used to make GaAs quan- 
tum well High Electron Mobility Transistors (HEMTs) which are used in all 
cellular telephones and in radio telescope receivers where they are prized for 
their low noise and ability to amplify extremely weak signals. The same tech- 
nology is widely utilized to produce the quantum well lasers used in compact 
disk players. 

1.1.3 Why is Disorder and Localization Important? 

Paradoxically, the extreme universality of the transport properties in the quan- 
tum Hall regime occurs because of, rather than in spite of, the random disorder 
and uncontrolled imperfections which the devices contain. Anderson localiza- 
tion in the presence of disorder plays an essential role in the quantization, but 
this localization is strongly modified by the strong magnetic field. 

In two dimensions (for zero magnetic field and non-interacting electrons) all 
states are localized even for arbitrarily weak disorder. The essence of this weak 
localization effect is the current 'echo' associated with the quantum interference 
corrections to classical transport [14]. These quantum interference effects rely 
crucially on the existence of time-reversal symmetry. In the presence of a strong 
quantizing magnetic field, time-reversal symmetry is destroyed and the localiza- 
tion properties of the disordered 2D electron gas are radically altered. We will 
shortly see that there exists a novel phase transition, not between a metal and 
insulator, but rather between two distinctly different insulating states. 

In the absence of any impurities the 2DEG is translationally invariant and 
there is no preferred frame of reference.^ As a result we can transform to a 
frame of reference moving with velocity —v relative to the lab frame. In this 
frame the electrons appear to be moving at velocity +v and carrying current 
density 

J = -nev, (1.4) 

where n is the areal density and we use the convention that the electron charge 
is — e. In the lab frame, the electromagnetic fields are 

^ = (1.5) 

■^This assumes that we can ignore the periodic potential of the crystal which is of course 
fixed in the lab frame. Within the efi'ective mass approximation this potential modifies the 
mass but does not destroy the Galilean invariance since the energy is still quadratic in the 
momentum. 



The Quantum Hall Effect 5 

B = Bz. (1.6) 
In the moving frame they are (to lowest order in v/c) 

E = --vxB (1.7) 
c 

B = Bz. (1.8) 



This Lorentz transformation picture is precisely equivalent to the usual state- 
ment that an electric field must exist which just cancels the Lorentz force ^vxB 
in order for the device to carry the current straight through without deflection. 
Thus we have 

E= —fx B 

nec 

The resistivity tensor is defined by 

Hence we can make the identification 

^ _B_ / +1\ 

= nec \ ~1 / 

The conductivity tensor is the matrix inverse of this so that 
and 

nec f ~1 \ 

Notice that, paradoxically, the system looks insulating since CTxx = and yet it 
looks like a perfect conductor since p^^ = 0. In an ordinary insulator a^y = 
and so p^x = oo. Here a^y = ^ and so the inverse exists. 

The argument given above relies only on Lorentz covariance. The only prop- 
erty of the 2DEG that entered was the density. The argument works equally well 
whether the system is classical or quantum, whether the electron state is liquid, 
vapor, or solid. It simply does not matter. Thus, in the absence of disorder, 
the Hall effect teaches us nothing about the system other than its density. The 
Hall resistivity is simply a linear function of magnetic field whose slope tells us 
about the density. In the quantum Hall regime we would therefore see none of 
the novel physics in the absence of disorder since disorder is needed to destroy 
translation invariance. Once the translation invariance is destroyed there is a 
preferred frame of reference and the Lorentz covariance argument given above 
fails. 

Figure (1.2) shows the remarkable transport data for a real device in the 
quantum Hall regime. Instead of a Hall resistivity which is simply a linear 



(1.9) 
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MAGNETIC FIELD (Tesia) 

Figure 1.2: Integer and fractional quantum Hall transport data showing the 
plateau regions in the Hall resistance i?H and associated dips in the dissipative 
resistance R. The numbers indicate the Landau level filling factors at which 
various features occur. After ref. [15]. 
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Figure 1.3: Persistent current circulating in a quantum Hall device having the 
Corbino geometry. The radial electric field is maintained by the charges which 
can not flow back together because a^x is nearly zero. Tliese cliargcs result 
from the radial current pulse associated with the azimuthal electric field pulse 
produced by the applied flux ^{t). 



function of magnetic field, we see a series of so-called Hall plateaus in which p^y 
is a universal constant 

P.y = -~ (1.14) 

independent of all microscopic details (including the precise value of the mag- 
netic field). Associated with each of these plateaus is a dramatic decrease in the 
dissipative resistivity pxx — which drops as much as 13 orders of magnitude 
in the plateau regions. Clearly the system is undergoing some sort of sequence 
of phase transitions into highly ideafized dissipationless states. Just as in a su- 
perconductor, the dissipationless state supports persistent currents. These can 
be produced in devices having the Corbino ring geometry shown in fig. (1.3). 
Applying additional fiux through the ring produces a temporary azimuthal elec- 
tric field by Faraday induction. A current pulse is induced at right angles to the 
E field and produces a radial charge polarization as shown. This polarization 
induces a (quasi-) permanent radial electric field which in turn causes persistent 
azimuthal currents. Torque magnetometer measurements [16] have shown that 
the currents can persist ^ 10^ sees at very low temperatures. After this time 
the tiny cr^x gradually allows the radial charge polarization to dissipate. We 
can think of the azimuthal currents as gradually spiraling outwards due to the 
Hall angle (between current and electric field) being very slightly less than 90° 
(by ~ 10-"). 

We have shown that the random impurity potential (and by implication 
Anderson localization) is a necessary condition for Hall plateaus to occur, but 
we have not yet understood precisely how this novel behavior comes about. 
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That is our next task. 



1.2 Classical and Semi- Classical Dynamics 
1.2.1 Classical Approximation 

The classical equations of motion for an electron of charge — e moving in two di- 
mensions rmdcr the influence of the Lorentz force ^vx B caused by a magnetic 
field B = Bz are 

eB 

mx = y (1-15) 

eB 

my = H X. (1-16) 

c 

The general solution of these equations corresponds to motion in a circle of 
arbitrary radius R 

f=R {coii{ujct + 6), sin(wci + S)) . (1.17) 

Here S is an arbitrary phase for the motion and 

eB 

^c=— 1-18 
mc 

is known as the classical cyclotron frequency. Notice that the period of the orbit 
is independent of the radius and that the tangential speed 

V = RjJc (1-19) 

controls the radius. A fast particle travels in a large circle but returns to the 
starting point in the same length of time as a slow particle which (necessarily) 
travels in a small circle. The motion is thus isochronous much like that of a 
harmonic oscillator whose period is independent of the amplitude of the motion. 
This apparent analogy is not an accident as we shall see when we study the 
Hamiltonian (which we will need for the full quantum solution). 

Because of some subtleties involving distinctions between canonical and me- 
chanical momentum in the presence of a magnetic field, it is worth reviewing 
the formal Lagrangian and Hamiltonian approaches to this problem. The above 
classical equations of motion follow from the Lagrangian 

£ = -mx^i^ - -x^A^', (1.20) 
2 c ^ ' 

where /i = 1,2 refers to x and y respectively and A is the vector potential 
evaluated at the position of the particle. (We use the Einstein summation 
convention throughout this discussion.) Using 

|±- = -V5,^/^ (1.21) 
ox" c 
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and 

5C e 

— = ma;" - -A"" (1.22) 
the Euler-Lagrange equation of motion becomes 

mx" = [d^A^" - a^A"] x^. (1.23) 

Using 

B = V X A (1.24) 

B°' = e"'f^'^dfiA'< (1.25) 

shows that this is equivalent to eqs. (1.15-1.16). 

Once we have the Lagrangian we can deduce the canonical momentum 

= mx''--A>', (1.26) 

c 

and the Hamiltonian 

H\p, x\ = x^^p" - £(:?, x) 

(Recall that the Lagrangian is canonically a function of the positions and ve- 
locities while the Hamiltonian is canonically a function of the positions and 
momenta). The quantity 

pLc^^P^ + I^^ (1-28) 
is known as the mechanical momentum. Hamilton's equations of motion 



p» = = -— (p- + i A") a^A" (1.30) 

ox^^ mc \ c I 



show that it is the mechanical momentum, not the canonical momentum, which 
is equal to the usual expression related to the velocity 

<ech = ^*''- (1-31) 

Using Hamilton's equations of motion we can recover Newton's law for the 
Lorentz force given in eq. (1-23) by simply taking a time derivative of x^ in 
eq. (1.29) and then using eq. (1.30). 

The distinction between canonical and mechanical momentum can lead to 
confusion. For example it is possible for the particle to have a finite velocity 
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while having zero (canonical) momentum! Furthermore the canonical momen- 
tum is dependent (as we will see later) on the choice of gauge for the vector 
potential and hence is not a physical observable. The mechanical momentum, 
being simply related to the velocity (and hence the current) is physically observ- 
able and gauge invariant. The classical equations of motion only involve the curl 
of the vector potential and so the particular gauge choice is not very important 
at the classical level. We will therefore delay discussion of gauge choices until 
we study the full quantum solution, where the issue is unavoidable. 



1.2.2 Semi-classical Approximation 

Recall that in the semi-classical approximation used in transport theory we 
consider wave packets j^^^-^ ^^^^ (r, t) made up of a linear superposition of Bloch 
waves. These packets are large on the; scale; of the de Broglie wavelength so that 
they have a well-defined central wave vector K{t), but they are small on the 
scale of everything else (external potentials, etc.) so that they simultaneously 
can be considered to have well-defined mean position R{t). (Note that K and R 
are parameters labeling the wave packet not arguments.) We then argue (and 
will discuss further below) that the solution of the Schrodinger equation in this 
semiclassical limit gives a wave packet whose parameters K{t) and R{t) obey 
the appropriate analog of the classical Hamilton equations of motion 

= ^fj, ^' (1-32) 
UK" = ' • (1.33) 

Naturally this leads to the same circular motion of the wave packet at the 
classical cyclotron frequency discussed above. For weak fields and fast electrons 
the radius of these circular orbits will be large compared to the size of the wave 
packets and the semi-classical approximation will be valid. However at strong 
fields, the approximation begins to break down because the orbits are too small 
and because becomes a significant (large) energy. Thus we anticipate that 
the semi-classical regime requires -C ep, where ep is the Fermi energy. 

We have already seen hints that the problem we are studying is really a 
harmonic oscillator problem. For the harmonic oscillator there is a characteristic 
energy scale Uu: (in this case ?iWc) and a characteristic length scale I for the zero- 
point fluctuations of the position in the ground state. The analog quantity in 
this problem is the so-called magnetic length 



(1.34) 

Itesla 




The physical interpretation of this length is that the area 27r^^ contains one 
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quantum of magnetic flux $o where^ 

$0 = (1.35) 
That is to say, the density of magnetic flux is 

B=||. (1.36) 

To be in the semiclassical hmit then requires that the Fermi wavelength be 
small on the scale of the magnetic length so that kpi '> 1. This condition turns 
out to be equivalent to hwc <C ep so they are not separate constraints. 



Exercise 1.1 Use the Bohr-Sommerfeld quantization condition that the orbit 

have a circumference containing an integral number of de Broglie wavelengths 
to find the allowed orbits of a 2D electron moving in a uniform magnetic field. 
Show that each successive orbit encloses precisely one additional quantum of 
flux in its interior. Hint: It is important to make the distinction between the 
canonical momentum (which controls the de Broglie wavelength) and the me- 
chanical momentum (which controls the velocity). The calculation is simplified 
if one uses the symmetric gauge A ~ ^\r x B in which the vector potential 
is purely azimuthal and independent of the azimuthal angle. 



1.3 Quantum Dynamics in Strong B Fields 

Since we will be dealing with the Hamiltonian and the Schrodinger equation, 
our first order of business is to choose a gauge for the vector potential. One 
convenient choice is the so-called Landau gauge: 

A{r) = xBy (1.37) 

which obeys V x ^ = Bz. In this gauge the vector potential points in the y 
direction but varies only with the x position, as illustrated in flg. (1.4). Hence 
the system still has translation invariance in the y direction. Notice that the 
magnetic field (and hence all the physics) is translationally invariant, but the 
Hamiltonian is not! (See exercise 1.2.) This is one of many peculiarities of 
dealing with vector potentials. 

Exercise 1.2 Show for the Landau gauge that even though the Hamiltonian 
is not invariant for translations in the x direction, the physics is still invariant 

since the change in the Hamiltonian that occurs under translation is simply 
equivalent to a gauge change. Prove this for any arbitrary gauge, assuming 
only that the magnetic field is uniform. 

^Note that in the study of superconductors the flux quantum is deflned with a factor of 2e 
rather than e to account for the pairing of the electrons in the condensate. 
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Figure 1.4: Illustration of the Landau gauge vector potential A = xBj). The 
magnetic field is perfectly uniform, but the vector potential has a preferred 
origin and orientation corresponding to the particular gauge choice. 



The Hamiltonian can be written in the Landau gauge as 

Taking advantage of the translation symmetry in the y direction, let us attempt 
a separation of variables by writing the wave function in the form 

Mx,y) = e'^y fk{x). (1.39) 

This has the advantage that it is an eigenstate of Py and hence we can make the 
replacement py — > Tik in the Hamiltonian. After separating variables we have 
the effective one-dimensional Schrodinger equation 

hkfk{x) = ekfk{x), (1.40) 

where 

1 2 1 . eB " 



This is simply a one-dimensional displaced harmonic oscillator^ 

hk = ^pl + \mu:l{x + k(:'f (1.42) 



^Thus we have arrived at the harmonic oscillator hinted at semiclassically, but paradoxically 
it is only one-dimensional, not two. The other degree of freedom appears (in this gauge) in 
the y momentum. 
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whose frequency is the classical cyclotron frcqiiency and whose central position 
= —k£^ is (somewhat paradoxically) determined by the y momentum quan- 
tum number. Thus for each plane wave chosen for the y direction there will be 
an entire family of energy eigenvalues 

ekn = {n+^)hwc (1.43) 

which depend only on n are completely independent of the y momentum hk. 
The corresponding (unnormalized) eigenfunctions are 

Mr) = -^e'^^H^ix + H2)e-ik(-+^^')^ (1.44) 
VL 

where Hn is (as usual for harmonic oscillators) the nth Hermite polynomial (in 
this case displaced to the now central position Xk). 



Exercise 1.3 Verify that eq. (1-44) is in fact a solution of the Schrddinger 
equation as claimed. 



These harmonic oscillator levels are called Landau levels. Due to the lack 
of dependence of the energy on k, the degeneracy of each level is enormous, as 
we will now show. We assume periodic boundary conditions in the y direction. 
Because of the vector potential, it is impossible to simultaneously have periodic 
boundary conditions in the x direction. However since the basis wave functions 
are harmonic oscillator polynomials multiplied by strongly converging gaussians, 
they rapidly vanish for positions away from the center position Xq = —kP. Let 
us suppose that the sample is rectangular with dimensions , Ly and that the 
left hand edge is at a: = — and the right hand edge is at a: = 0. Then the 
values of the wavevector k for which the basis state is substantially inside the 
sample run from fc = to fc = L^jH?' . It is clear that the states at the left edge 
and the right edge differ strongly in their k values and hence periodic boundary 
conditions are impossible.^ 

The total number of states in each Landau level is then 

where 

^ ^ (1.46) 

is the number of flux quanta penetrating the sample. Thus there is one state per 
Landau level per flux quantum which is consistent with the semiclassical result 
from Exercise (1.1). Notice that even though the family of allowed wavevectors 

^The best one can achieve is so-called quasi-periodic boundary conditions in which the 
phase difference between the left and right edges is zero at the bottom and rises linearly 
with height, reaching = LxLy/£^ at the top. The eigenfunctions with these boundary 

conditions are elliptic theta functions which are linear combinations of the gaussians discussed 
here. See the discussion by Haldane in Ref. [3]. 
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is only onc-dimcnsional, wc find that the degeneracy of each Landau level is 
extensive in the two-dimensional area. The reason for this is that the spacing 
between wave vectors allowed by the periodic boundary conditions = |^ 
decreases while the range of allowed wave vectors [0, Lx/t"^] increases with in- 
creasing L. The reader may also worry that for very large samples, the range 
of allowed values of k will be so large that it will fall outside the first Brillouin 
zone forcing us to inchidc band mixing and the periodic lattice potential beyond 
the effective mass approximation. This is not true however, since the canonical 
momentum is a gauge dependent quantity. The value of k in any particular 
region of the sample can be made small by shifting the origin of the coordinate 
system to that region (thereby making a gauge transformation). 

The width of the harmonic oscillator wave functions in the nth Landau level 
is of order ^fnl. This is microscopic compared to the system size, but note that 
the spacing between the centers 

A = Aa' = ^ (1.47) 

is vastly smaller (assuming Ly » £). Thus the supports of the different basis 

states are strongly o\^erlappiiig (l)ut tlixn- nvv still ortliogoual). 



Exercise 1.4 Using the fact that the energy for the nth Juirmonie oscillator 
state is {n + ^)?iUc, present a semi-classical argument explaining the result 
claimed above that the width of the support of the wave function scales as 

Exercise 1.5 Using the Landau gauge, construct a gaussian wave packet in 
the lowest Landau level of the form 

+00 



*(a;,i/) = / akc'^ye 

J —oo 

choosing ak in such a way that the wave packet is localized as closely as possible 
around some point R. What is the smallest size wave packet that can be 
constructed without mixing in higher Landau levels? 



Having now found the cigcnfunctions for an electron in a strong magnetic 
field we can relate them back to the semi-classical picture of wave packets under- 
going circular cyclotron motion. Consider an initial semiclassical wave packet 
located at some position and having some specified momentum. In the semi- 
classical limit the mean energy of this packet will greatly exceed the cyclotron 
energy 3> hujc and hence it will be made up of a linear combination of a 

large number of different Landau level states centered around n — ^ ^ 



nr,t) =Y,j a„(^)V„fe(r )e-^("+^)-^*. (1.48) 

Notice that in an ordinary 2D problem at zero field, the complete set of plane 
wave states would be labeled by a 2D continuous momentum label. Here we 



The Quantum Hall Effect 



15 



have one discrete label (the Landau level index) and a ID continuous labels 
(the y wave vector). Thus the 'sum' over the complete set of states is actually 
a combination of a summation and an integration. 

The details of the initial position and momentum are controlled by the am- 
plitudes an{k). We can immediately see however, that since the energy levels 
are exactly evenly spaced that the motion is exactly periodic: 

*(r,f+ — ) = *(r,t). (1.49) 

If one works through the details, one finds that the motion is indeed circular 
and corresponds to the expected semi-classical cyclotron orbit. 

For simplicity we will restrict the remainder of our discussion to the lowest 
Landau level where the (correctly normalized) cigcnfunctions in the Landau 
gauge are (dropping the index n = from now on): 

V'fe(f ) = ^^^e'^ye-M^+^^')' (1.50) 

and every state has the same energy eigenvalue = ^tiujc- 

We imagine that the magnetic field (and hence the Landau level splitting) is 
very large so that we can ignore higher Landau levels. (There are some subtleties 
here to which we will return.) Because the states are all degenerate, any wave 
packet made up of any combination of the basis states will be a stationary 
state. The total current will therefore be zero. We anticipate however from 
semiclassical considerations that there should be some remnant of the classical 
circular motion visible in the local current density. To see this note that the 
expectation value of the current in the fcth basis state is 

(/) = -e- hk (p+ -A) ^k) . (1.51) 
m \ \ c / / 



The y component of the current is 



eB \ — 1 (2,+fe^2)2 
nk H X \ e ' 



= + (1-52) 

We see from the integrand that the current density is antisymmetric about the 
peak of the gaussian and hence the total current vanishes. This antisymmetry 
(positive vertical current on the left, negative vertical current on the right) is 
the remnant of the semiclassical circular motion. 

Let us now consider the case of a uniform electric field pointing in the x 
direction and giving rise to the potential energy 

V{r) = +eEx. (1.53) 



This still has translation symmetry in the y direction and so our Landau gauge 
choice is still the most convenient. Again separating variables we see that the 
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solution is nearly the same as before, except that the displacement of the har- 
monic oscillator is slightly different. The Hamiltonian in eq. (1-54) becomes 

hk = + ^mujl {x + kff + eEx. (1.54) 

Completing the square we see that the oscillator is now centered at the new 
position 

Xk = -ke - (1.55) 

and the energy eigenvalue is now linearly dependent on the particle's peak po- 
sition Xk (and therefore linear in the y momentum) 

e/c = ^?iwc + e-EXfc + imw^, (1.56) 

where 

V = -c|. (1.57) 

Because of the shift in the peak position of the wavefunction, the perfect anti- 
symmetry of the current distribution is destroyed and there is a net current 

{Jy) = -ev (1.58) 

showing that vy is simply the usual cE x B /B^ drift velocity. This result can 
be derived either by explicitly doing the integral for the current or by noting 
that the wave packet group velocity is 

1 dek eE dXk _ 

independent of the value of k (since the electric field is a constant in this case, 
giving rise to a strictly linear potential). Thus we have recovered the correct 
kinematics from our quantum solution. 

It should be noted that the applied electric field 'tilts' the Landau levels in 
the sense that their energy is now linear in position as illustrated in fig. (1.5). 
This means that there are degeneracies between different Landau level states 
because different kinetic energy can compensate different potential energy in 
the electric field. Nevertheless, we have found the exact eigenstates (i.e., the 
stationary states). It is not possible for an electron to decay into one of the other 
degenerate states because they have different canonical momenta. If however 
disorder or phonons are available to break translation symmetry, then these 
decays become allowed and dissipation can appear. The matrix elements for 
such processes are small if the electric field is weak because the degenerate 
states are widely separated spatially due to the small tilt of the Landau levels. 
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(a) (b) 

Figure 1.5: Illustration of electron Landau energy levels (n + i) TicUc vs. position 
Xk = ~k£'^- (a) Zero electric field case, (b) Case with finite electric field pointing 
in the +x direction. 



Exercise 1.6 It is interesting to note that the exact eigenstates in the pres- 
ence of the electric field can be viewed as displaced oscillator states in the 
original (zero E field) basis. In this basis the displaced states are linear com- 
binations of all the Landau level excited states of the sarn,e k. Use first-order 
perturbation theory to find the amount by which the n = 1 Landau level is 
mixed into the n = state. Compare this with the exact amount of mixing 
computed using the exact displaced oscillator state. Show that the two results 
agree to first order in E. Because the displaced state is a linear combination of 
more than one Landau level, it can carry a finite current. Give an argument, 
based on perturbation theory why the amount of this current is inversely pro- 
portional to the B field, but is independent of the mass of the particle. Hint: 
how does the mass affect the Landau level energy spacing and the current 
operator? 



1.4 IQHE Edge States 

Now that we understand drift in a uniform electric field, we can consider the 

problem of electrons confined in a Hall bar of finite width by a non-uniform 
electric field. For simplicity, we will consider the situation where the potential 
V{x) is smooth on the scale of the magnetic length, but this is not central to 
the discussion. If we assume that the system still has translation symmetry in 
the y direction, the solution to the Schrodinger equation must still be of the 
form 



Jp{x,y) = 




(1.60) 



The function fk will no longer be a simple harmonic wave function as we found 
in the case of the uniform electric field. However we can anticipate that will 
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Figure 1.6: Illustration of a smooth confining potential which varies only in the 
X direction. The horizontal dashed line indicates the equilibrium fermi level. 

The dashed curve indicates the wave packet envelope fk which is displaced from 
its nominal position Xk = —ki^ by the slope of the potential. 



still be peaked near (but in general not precisely at) the point = —ki"^. The 
eigenvalues will no longer be precisely linear in k but will still reflect the 
kinetic energy of the cyclotron motion plus the local potential energy V{Xk) 
(plus small corrections analogous to the one in eq. (1.56)). This is illustrated in 
fig. (1.6). We see that the group velocity 

"^ = nafc^ ^'-''^ 

has the opposite sign on the two edges of the sample. This means that in the 
ground state there are edge currents of opposite sign flowing in the sample. The 
semi-classical interpretation of these currents is that they represent 'skipping 
orbits' in which the circular cyclotron motion is interrupted by collisions with 
the walls at the edges as illustrated in flg. (1.7). 

One way to analyze the Hall effect in this system is quite analogous to the 
Landauer picture of transport in narrow wires [17, 18]. The edge states play 
the role of the left and right moving states at the two fermi points. Because 
(as we saw earlier) momentum in a magnetic field corresponds to position, the 
edge states are essentially real space realizations of the fermi surface. A Hall 
voltage drop across the sample in the x direction corresponds to a difference in 
electrochemical potential between the two edges. Borrowing from the Landauer 
formulation of transport, we will choose to apply this in the form of a chemical 
potential difference and ignore any changes in electrostatic potential.^ What 
this does is increase the number of electrons in skipping orbits on one edge of 
the sample and/or decrease the number on the other edge. Previously the net 
current due to the two edges was zero, but now there is a net Hall current. To 
calculate this current we have to add up the group velocities of all the occupied 



^This has led to various confusions in the hterature. If there is an electrostatic potential 
gradient then some of the net Hall current may be carried in the bulk rather than at the 
edges, but the final answer is the same. In any case, the essential part of the physics is that 
the only place where there are low lying excitations is at the edges. 
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Figure 1.7: Semi-classical view of skipping orbits at the fermi level at the two 
edges of the sample where the confining electric field causes E x B drift. The 
circular orbit illustrated in the center of the sample carries no net drift current 
if the local electric field is zero. 



states 



where for the moment we assume that in the bulk, only a single Landau level is 

occupied and is the probability that state k in that Landau level is occupied. 
Assuming zero temperature and noting that the integrand is a perfect derivative, 
we have 

^=-^''''c!e = -|[ML-/Xfl]. (1.63) 

(To understand the order of limits of integration, recall that as k increases, Xk 
decreases.) The definition of the Hall voltage drop is"" 

{+e)VH = (+e) [Vr - Vl] = [i^r - hl] ■ (1-64) 

Hence 

/ = -^jVh, (1.65) 

where we have now allowed for the possibility that u different Landau levels are 
occupied in the bulk and hence there are v separate edge channels contributing 

'^To get the signs straight here, note that an increase in chemical potential brings in more 
electrons. This is equivalent to a more positive voltage and hence a more negative potential 
energy —eV. Since H — ij,N enters the thermodynamics, electrostatic potential energy and 
chemical potential move the electron density oppositely. V and // thus have the same sign of 
effect because electrons are negatively charged. 
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to the current. This is the analog of having i/ 'open' channels in the Landauer 
transport picture. In the Landauer picture for an ordinary wire, we are consid- 
ering the longitudinal voltage drop (and computing a^x), while here we have 
the Hall voltage drop (and arc computing (Jxy)- The analogy is quite precise 
however because we view the right and left movers as having distributions con- 
trolled by separate chemical potentials. It just happens in the QHE case that 
the right and left movers are physically separated in such a way that the voltage 
drop is transverse to the current. Using the above result and the fact that the 
current flows at right angles to the voltage drop we have the desired results 

(Txx = (1.66) 



'xy 



(1.67) 



with the quantum number i/ being an integer. 

So far wo have been ignoring the possible effects of disorder. Recall that 
for a single-channel one-dimensional wire in the Landauer picture, a disordered 
region in the middle of the wire will reduce the conductivity to 

I=j\T\^ (1.68) 

where |Tp is the probability for an electron to be transmitted through the 
disordered region. The reduction in transmitted current is due to back scattering. 
Remarkably, in the QHE case, the back scattering is essentially zero in very 
wide samples. To see this note that in the case of the Hall bar, scattering into 
a backward moving state would require transfer of the electron from one edge 
of the sample to the other since the edge states arc spatially separated. For 
samples which are very wide compared to the magnetic length (more precisely, 
to the Anderson localization length) the matrix element for this is exponentially 
small. In short, there can be nothing but forward scattering. An incoming wave 
given by eq. (1.60) can only be transmitted in the forward direction, at most 
suffering a simple phase shift 6k 

Vout(a;,y) = -^e^'^e'^yhix). (1.69) 

This is because no other states of the same energy are available. If the disorder 
causes Landau level mixing at the edges to occur (because the confining potential 
is relatively steep) then it is possible for an electron in one edge channel to 
scatter into another, but the current is still going in the same direction so 
that there is no reduction in overall transmission probability. It is this chiral 
(unidirectional) nature of the edge states which is responsible for the fact that 
the Hall conductance is correctly quantized independent of the disorder. 

Disorder will broaden the Landau levels in the bulk and provide a reservoir 
of (localized) states which will allow the chemical potential to vary smoothly 
with density. These localized states will not contribute to the transport and so 
the Hall conductance will be quantized over a plateau of finite width in B (or 
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density) as seen in the data. Thus obtaining the universal vahie of quantized 
Hall conductance to a precision of 10^^° does not require fine tuning the applied 
B field to a similar precision. 

The localization of states in the bulk by disorder is an essential part of 
the physics of the quantum Hall effect as we saw when we studied the role of 
translation invariance. We learned previously that in zero magnetic field all 
states are (weakly) localized in two dimensions. In the presence of a quantizing 
magnetic field, most states are strongly localized as discussed above. However 
if all states were localized then it would be impossible to have a quantum phase 
transition from one QHE plateau to the next. To understand how this works 
it is convenient to work in a semiclassical percolation picture to be described 
below. 

Exercise 1.7 Show that the number of edge channels whose energies lie in 
the gap between two Landau levels scales with the length L of the sample, while 
the number of bulk states scales with the area. Use these facts to show that 
the range of magnetic field in which the chemical potential lies in between two 
Landau levels scales to zero in the thermodynamic limit. Hence finite width 
quantized Hall plateaus can not occur in the absence of disorder that produces 
a reservoir of localized states in the bulk whose number is proportional to the 
area. 



1.5 Semiclassical Percolation Picture 

Let us consider a smooth random potential caused, say, by ionized silicon donors 
remotely located away from the 2DEG in the GaAs semiconductor host. We 
take the magnetic field to be very large so that the magnetic length is small on 
the scale over which the potential varies. In addition, we ignore the Coulomb 
interactions among the electrons. 

What is the nature of the eigenfunctions in this random potential? We have 
learned how to solve the problem exactly for the case of a constant electric field 
and know the general form of the solution when there is translation invariance in 
one direction. We found that the wave functions were plane waves running along 
lines of constant potential energy and having a width perpendicular to this which 
is very small and on the order of the magnetic length. The reason for this is the 
discreteness of the kinetic energy in a strong magnetic field. It is impossible for 
an electron stuck in a given Landau level to continuously vary its kinetic energy. 
Hence energy conservation restricts its motion to regions of constant potential 
energy. In the limit of infinite magnetic field where Landau level mixing is 
completely negligible, this confinement to lines of constant potential becomes 
exact (as the magnetic length goes to zero). 

We are led to the following somewhat paradoxical picture. The strong mag- 
netic field should be viewed as putting the system in the quantum limit in the 
sense that hUc is a very large energy (comparable to ep)- At the same time 
(if one assumes the potential is smooth) one can argue that since the magnetic 
length is small compared to the scale over which the random potential varies, 
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the system is in a semi-classical limit where small wave packets (on the scale of 
£) follow classical E x B drift trajectories. 

Prom this discussion it then seems very reasonable that in the presence 
of a smooth random potential, with no particular translation symmetry, the 
eigenfunctions will live on contour lines of constant energy on the random energy 
surface. Thus low energy states will be found lying along contours in deep valleys 
in the potential landscape while high energy states will be found encircling 
'mountain tops' in the landscape. Naturally these extreme states will be strongly 
localized about these extrema in the potential. 



Exercise 1.8 Using the Lagrangian for a charged particle in a magnetic field 
with a scalar potential V{f ), consider the high field limit by setting the mass 
to zero (thereby sending the quantum cyclotron energy to infinity). 

1. Derive the classical equations of motion from the Lagrangian and show 
that they yield simple E x B drift along isopotential contours. 

2. Find the momentum conjugate to the coordinate x and show that (with 
an appropriate gauge choice) it is the coordinate y: 



In the infinite field limit where £ ^ the coordinates commute and we re- 
cover the semi-classical result in which effectively point particles drift along 
isopotentials. 



To understand the nature of states at intermediate energies, it is useful 
to imagine gradually filling a random landscape with water as illustrated in 
fig. (1.8). In this analogy, sea level represents the chemical potential for the 
electrons. When only a small amount of water has been added, the water will 
fill the deepest valleys and form small lakes. As the sea level is increased the 
lakes will grow larger and their shorelines will begin to take on more complex 
shapes. At a certain critical value of sea level a phase transition will occur in 
which the shoreline percolates from one side of the system to the other. As the 
sea level is raised still further, the ocean will cover the majority of the land and 
only a few mountain tops will stick out above the water. The shore line will no 
longer percolate but only surround the mountain tops. 

As the sea level is raised still higher additional percolation transitions will 
occur successively as each successive Landau level passes under water. If Lan- 
dau level mixing is small and the disord(;r potential is symmetrically distributed 
about zero, then the critical value of the chemical potential for the nth perco- 



n 

Px = --^y 



(1.70) 



so that we have the strange commutation relation 



[x,y] = -ie. 



(1.71) 
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Figure 1.8: Contour map of a smooth random landscape. Closed dashed lines 
indicate local mountain peaks. Closed solid lines indicate valleys. From top to 
bottom, the gray filled areas indicate the increasing 'sea level' whose shoreline 
finally percolates from one edge of the sample to the other (bottom panel) . The 
particle-hole excitations live along the shoreline and become gapless when the 
shoreline becomes infinite in extent. 

lation transition will occur near the center of the nth Landau level 

fi: = {n+^)nue. (1.72) 

This percolation transition corresponds to the transition between quantized 
Hall plateaus. To sec why, note that when the sea level is below tlic percolation 
point, most of the sample is dry land. The electron gas is therefore insulating. 
When sea level is above the percolation point, most of the sample is covered 
with water. The electron gas is therefore connected throughout the majority of 
the sample and a quantized Hall current can be carried. Another way to see this 
is to note that when the sea level is above the percolation point, the confining 
potential will make a shoreline along the full length of each edge of the sample. 
The edge states will then carry current from one end of the sample to the other. 

We can also understand from this picture why the dissipative conductivity 
a^x bas a sharp peak just as the plateau transition occurs. (Recall the data 
in fig. (1.2).) Away from the critical point the circumference of any particular 
patch of shoreline is finite. The period of the semiclassical orbit around this is 
finite and hence so is the quantum level spacing. Thus there are small energy 
gaps for excitation of states across these real-space fermi levels. Adding an 
infinitesimal electric field will only weakly perturb these states due to the gap 
and the finiteness of the perturbing matrix element which will be limited to 
values on the order of ~ eED where D is the diameter of the orbit. If however 
the shoreline percolates from one end of the sample to the other then the orbital 
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Figure 1.9: Illustration of edge states that wander deep into the bulk as the 
quantum Hall localization transition is approached from the conducting side. 
Solid arrows indicate the direction of drift along the isopotential lines. Dashed 
arrows indicate quantum tunneling from one semi-classical orbit (edge state) to 
the other. This backscattering localizes the eigenstates and prevents transmis- 
sion through the sample using the 'edge' states (which become part of the bulk 
localized states). 



period diverges and the gap vanishes. An infinitesimal electric field can then 
cause dissipation of energy. 

Another way to see this is that as the percolation level is approached from 
above, the edge states on the two sides will begin taking detours deeper and 
deeper into the bulk and begin communicating with each other as the local- 
ization length diverges and the shoreline zig zags throughout the bulk of the 
sample. Thus electrons in one edge state can be back scattered into the other 
edge states and ultimately reflected from the sample as illustrated in fig. (1.9). 

Because the random potential broadens out the Landau level density of 
states, the quantized Hall plateaus will have finite width. As the chemical 
potential is varied in the regime of localized states in between the Landau level 
peaks, only the occupancy of localized states is changing. Hence the transport 
properties remain constant until the next percolation transition occurs. It is 
important to have the disorder present to produce this finite density of states 
and to localize those states. 

It is known that as the (classical) percolation point is approached in two 
dimensions, the characteristic size (diameter) of the shoreline orbits diverges 
like 

e ~ (1-73) 

where 5 measures the deviation of the sea level from its critical value. The 
shoreline structure is not smooth and in fact its circumference diverges with a 
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larger exponent 7/3 showing that these arc highly ramified fractal objects whose 
circumference scales as the 7/4th power of the diameter. 

So far we have assumed that the magnetic length is essentially zero. That is, 
we have ignored the fact that the wave function support extends a small distance 
transverse to the isopotential lines. If two different orbits with the same energy 
pass near each other but are classically disconnected, the particle can still tunnel 
between them if the magnetic length is finite. This quantum tunneling causes the 
localization length to diverge faster than the classical percolation model predicts. 
Numerical simulations find that the localization length diverges like [19-22] 



where the exponent ly (not to be confused with the Landau level filling factor!) 
has a value close (but probably not exactly equal to) 7/3 rather than the 4/3 
found in classical percolation. It is believed that this exponent is universal and 
independent of Landau level index. 

Experiments on the quantum critical behavior arc quite difficult but there 
is evidence [23], at least in selected samples which show good scaling, that ly is 
indeed close to 7/3 (although there is some recent controversy on this point. [24]) 
and that the conductivity tensor is universal at the critical point. [21,25] Why 
Coulomb interactions that are present in real samples do not spoil agreement 
with the numerical simulations is something of a mystery at the time of this 
writing. For a discussion of some of these issues see [13]. 



Under some circiimstanccs of weak (but non-zero) disorder, quantized Hall 
plateaus appear which arc characterized by simple rational fractional quantum 
numbers. For example, at magnetic fields three times larger than those at which 
the 1^=1 integer filling factor plateau occurs, the lowest Landau level is only 1/3 
occupied. The system ought to be below the percolation threshold and hence 
be insulating. Instead a robust quantized Hall plateau is observed indicating 
that electrons can travel through the sample and that (since axx — > 0) there 
is an excitation gap. This novel and quite unexpected physics is controlled by 
Coulomb repulsion between the electrons. It is best understood by first ignoring 
the disorder and trying to discover the nature of the special correlated many- 
body ground state into which the electrons condense when the filling factor is 
a rational fraction. 

For reasons that will become clear later, it is convenient to analyze the 
problem in a new gauge 



known as the symmetric gauge. Unlike the Landau gauge which preserves trans- 
lation symmetry in one direction, the symmetric gauge preserves rotational sym- 
metry about the origin. Hence we anticipate that angular momentum (rather 
than y linear momentum) will be a good quantum number in this gauge. 



(1.74) 



1.6 Fractional QHE 




(1.75) 
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For simplicity wc will restrict our attention to the lowest Landau level only 
and (simply to avoid some awkward minus signs) change the sign of the B field: 
B = —Bz. With these restrictions, it is not hard to show that the solutions of 
the free-particle Schrodinger equation having definite angular momentum are 

<fim = , ^ ^'"e-^l^l' (1.76) 

where z = (x+iy) /£ is a dimensionless complex number representing the position 

vector r = (x, y) and to > is an integer. 



Exercise 1.9 Verify that the basis functions in eq. (1.76) do solve the 
Schrodinger equation in the absence of a potential and do lie in the lowest 
Landau level. Hint: Rewrite the kinetic energy in such a way that p- A be- 
comes B ■ L. 



The angular momentum of these basis states is of course Urn. If we restrict 
our attention to the lowest Landau level, then there exists only one state with 
any given angular momentum and only non-negative values of to are allowed. 
This 'handedness' is a result of the chirality built into the problem by the 
magnetic field. 

It seems rather peculiar that in the Landau gauge we had a continuous one- 
dimensional family of basis states for this two-dimensional problem. Now we 
find that in a diff'erent gauge, we have a discrete one dimensional label for the 
basis states! Nevertheless, we still end up with the correct density of states 
per unit area. To see this note that the peak value of |<PmP occurs at a radius 
of i?peak = V2mP. The area 2-K(^m of a circle of this radius contains m flux 
quanta. Hence we obtain the standard result of one state per Landau level per 
quantum of flux penetrating the sample. 

Because all the basis states are degenerate, any linear combination of them 
is also an allowed solution of the Schrodinger equation. Hence any function of 
the form [26] 

*(a;,y) = /(^)e-^l^l' (1.77) 

is allowed so long as / is analytic in its argument. In particular, arbitrary 
polynomials of any degree N 

N 

f{z) = \l{z-Z,) (1.78) 

are allowed (at least in the thermodynamic limit) and are conveniently defined 
by the locations of their N zeros {Z^; j = 1, 2, . . . , iV}. 

Another useful solution is the so-called coherent state which is a particular 
infinite order polynomial 
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The wave function using this polynomial has the property that it is a narrow 
gaussian wave packet centered at the position defined by the complex number 
A. Completing the square shows that the probability density is given by 

|vl.,|2 = |/,|V^N^ = J^e-il-^l^ (1.80) 

This is the smallest wave packet that can be constructed from states within the 
lowest Landau level. The reader will find it instructive to compare this gaussian 
packet to the one constructed in the Landau gauge in exercise (1.5). 

Because the kinetic energy is completely degenerate, the effect of Coulomb 
interactions among the particles is nontrivial. To develop a feel for the problem, 
let us begin by solving the two-body problem. Recall that the standard proce- 
dure is to take advantage of the rotational symmetry to write down a solution 
with the relative angular momentum of the particles being a good quantum 
number and then solve the Schrodinger equation for the radial part of the wave 
function. Here we find that the analyticity properties of the wave functions in 
the lowest Landau level greatly simplifies the situation. If we know the angular 
behavior of a wave function, analyticity uniquely defines the radial behavior. 
Thus for example for a single particle, knowing that the angular part of the 
wave function is e""" , we know that the full wave function is guaranteed to 

o 1 I 1 2 1 I 1 2 

uniquely be r™e e~*'^' = z"^e~^'^' . 

Consider now the two body problem for particles with relative angular mo- 
mentum m and center of mass angular momentum M. The unique analytic 
wave function is (ignoring normalization factors) 

'bmM{Zl,Z2) = [Zi - 22)™(Z1 + z^)^ e--^^\'^\"+\'-\"\ (1.81) 

If m and M are no n- negative integers, then the prefactor of the exponential is 
simply a polynomial in the two arguments and so is a state made up of linear 
combinations of the degenerate one-body basis states given in eq. (1.76) 
and therefore lies in the lowest Landau level. Note that if the particles are 
spinless fermions then m must be odd to give the correct exchange symmetry. 
Remarkably, this is the exact (neglecting Landau level mixing) solution for the 
Schrodinger equation for any central potential y(|zi — Z2I) acting between the 
two particles.^ We do not need to solve any radial equation because of the 
powerful restrictions due to analyticity. There is only one state in the (lowest 
Landau level) Hilbert space with relative angular momentum m and center of 
mass angular momentum M . Hence (neglecting Landau level mixing) it is an 
exact eigenstate of any central potential, is the exact answer independent 

of the Hamiltonian! 

The corresponding energy eigenvalue Vm is independent of M and is referred 
to as the mth Haldane pseudopotential 

_ {mM\V\mM) 

{mM\mM) ' ^''^^^ 



^Note that neglecting Landau level mixing is a poor approximation for strong potentials 
V S> huic unless they are very smooth on the scale of the magnetic length. 
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Figure 1.10: The Haldane pseudopotcntial Vm vs. relative angidar momentum m 
for two particles interacting via the Coulomb interaction. Units are /e£, where 
e is the dielectric constant of the host semiconductor and the finite thickness of 
the quantum well has been neglected. 



The Haldane pseudopotentials for the repulsive Coulomb potential are shown 
in fig. (1.10). These discrete energy eigenstates represent bound states of the 
repulsive potential. If there were no magnetic field present, a repulsive potential 
would of course have only a continuous spectrum with no discrete bound states. 
However in the presence of the magnetic field, there are effectively bound states 
because the kinetic energy has been quenched. Ordinarily two particles that 
have a lot of potential energy because of their repulsive interaction can fly apart 
converting that potential energy into kinetic energy. Here however (neglecting 
Landau level mixing) the particles all have fixed kinetic energy. Hence particles 
that are repelling each other are stuck and can not escape from each other. 
One can view this semi-classically as the two particles orbiting each other under 
the influence of E x B drift with the Lorentz force preventing them from flying 
apart. In the presence of an attractive potential the eigenvalues change sign, but 
of course the eigenfunctions remain exactly the same (since they are unique)! 

The fact that a repulsive potential has a discrete spectrum for a pair of 
particles is (as we will shortly see) the central feature of the physics under- 
lying the existence of an excitation gap in the fractional quantum Hall effect. 
One might hope that since we have found analyticity to uniquely determine 
the two-body eigenstates, we might be able to determine many-particle eigen- 
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states exactly. The situation is complicated however by the fact that for three 
or more particles, the various relative angular momenta ^12,^13,^23, etc. do 
not all commute. Thus we can not write down general exact eigenstates. We 
will however be able to use the analyticity to great advantage and make exact 
statements for certain special cases. 



Exercise 1.10 Express the exact lowest Landau level two-body eigenstate 

in terms of the basis of all possible two-body Slater determinants. 

Exercise 1.11 Verify the claim that the Haldane pseudopotential Vm is inde- 
pendent of the center of mass angular momentum M. 

Exercise 1.12 Evaluate the Haldane pseudopotentials for the Coulomb po- 

2 2 
tential ^ . Express your answer in units of ^ . For the specific case 0/ e = 10 

and B = 10 T, express your answer in Kelvin. 

Exercise 1.13 Take into account the finite thickness of the quantum well by 
assuming that the one-particle basis states have the form 

where s is the coordinate in the direction normal to the quantum well. Write 
down (hut do not evaluate) the formal expression for the Haldane pseudo- 
potentials in this case. Qualitatively describe the effect of finite thickness on 
the values of the different pseudopotentials for the case where the well thickness 
is approximately equal to the magnetic length. 



1.6.1 The v —1 many-body state 

So far we have found the one- and two-body states. Our next task is to write 
down the wave function for a fully filled Landau level. We need to find 



m = m 



(1.83) 



where [z] stands for (zi, Z2, . . . , zjv) and / is a polynomial representing the 
Slater determinant with all states occupied. Consider the simple example of two 
particles. We want one particle in the orbital <^o and one in t^i, as illustrated 
schematically in fig. (1.11a). Thus (again ignoring normalization) 



m = 



{Z2 - Zi) 



= {zif{Z2f-{z2)\z,f 



(1.84) 
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m=0 m=l m=2 m=3 m=4 

Figure 1.11: Orbital occupancies for the maximal density filled Landau level 
state with (a) two particles and (b) three particles. There are no particle labels 
here. In the Slater determinant wave function, the particles are labeled but a 
sum is taken over all possible permutations of the labels in order to antisym- 
metrize the wave function. 



This is the lowest possible order polynomial that is antisymmetric. For the case 
of three particles we have (see fig. (1.11b)) 



{Zl)° {Z2f {Z,f 

{z^Y {z^y {zsY 

(zi)' {Z2f (zs)' 

-{Zl - Z2){zi - Z3){Z2 
3 

■YliZi-Zj) 

i<j 



Z2ZI 



Z-iZi 



'1^3 



ZzZ{ + 



Z\Z2 ■ 



'2^1 



(1.85) 



This form for the Slater determinant is known as the Vandermonde polynomial. 

The overall minus sign is unimportant and we will drop it. 

The single Slater determinant to fill the first N angular momentum states 
is a simple generalization of eq. (1.85) 



JV 



InM = Y[{zi - zj). 



(1.86) 



To prove that this is true for general N, note that the polynomial is fully 
antisymmetric and the highest power of any z that appears is Thus the 

highest angular momentum state that is occupied is m = A'' — 1. But since the 
antisymmetry guarantees that no two particles can be in the same state, all N 
states from m = Otom = A''— 1 must be occupied. This proves that we have 
the correct Slater determinant. 



Exercise 1.14 Show carefully that the Vandermonde polynomial for N par- 
ticles is in fact totally antisymmetric. 



One can also use induction to show that the Vandermonde polynomial is the 
correct Slater determinant by writing 



JV 



fN+l{z) = fN{z) YliZi - ZN+l) 



(1.87) 



The Quantum Hall Effect 



31 



which can be shown to agree with the resuh of expanding the determinant of 
the (A^+1) X matrix in terms of the minors associated with the (-/V + l)st 

row or column. 

Note that since the Vandermonde polynomial corresponds to the filled Lan- 
dau level it is the unique state having the maximum density and hence is an 
exact eigenstate for any form of interaction among the particles (neglecting 
Landau level mixing and ignoring the degeneracy in the center of mass angular 
momentum) . 

The (unnormalized) probability distribution for particles in the filled Landau 
level state is 

N 

= n l^i - ^jf e"*^"=i (1.88) 

i<j 

This seems like a rather complicated object about which it is hard to make any 
useful statements. It is clear that the polynomial term tries to keep the particles 
away from each other and gets larger as the particles spread out. It is also clear 
that the exponential term is small if the particles spread out too much. Such 
simple questions as, 'Is the density uniform?', seem hard to answer however. 

It turns out that there is a beautiful analogy to plasma physics developed 
by R. B. Laughlin which sheds a great deal of light on the nature of this many 
particle probability distribution. To see how this works, let us pretend that the 
norm of the wave function 



Z = J d'zi...J d'zN IVwr (1-89) 

is the partition function of a classical statistical mechanics problem with Boltz- 
mann weight 

\^[z]\^ = e-^^'^'-"" (1.90) 
where = — and 

[/class = m^Y.{~ \n\zi - Zj\) + ^Y1 l^'^l^- (1-91) 

i<j k 

(The parameter to = 1 in the present case but we introduce it for later conve- 
nience.) It is perhaps not obvious at first glance that we have made tremendous 
progress, but we have. This is because C/ciass turns out to be the potential energy 
of a fake classical one-component plasma of particles of charge to. in a miiform 
('jellium') neutralizing background. Hence we can bring to bear well-developed 
intuition about classical plasma physics to study the properties of l^'p. Please 
remember however that all the statements we make here are about a partic- 
ular wave function. There are no actual long-range logarithmic interactions 
in the quantum Hamiltonian for which this wave function is the approximate 
groundstate. 
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To imdcrstand this, let us first review the electrostatics of charges in three 
dimensions. For a charge Q particle in 3D, the surface integral of the electric 
field on a sphere of radius R surrounding the charge obeys 



dA-E = 4TrQ. (1.92) 

Since the area of the sphere is AttR'^ we deduce 

E{r) = (1.93) 

ip(f) = — (1.94) 
r 

and 

V • ^ = -VV = 47rQ 6^{f) (1.95) 

where <f is the electrostatic potential. Now consider a two-dimensional world 

where all the field lines arc confined to a plane (or equivalcntly consider the 
electrostatics of infinitely long charged rods in 3D). The analogous equation for 
the line integral of the normal electric field on a circle of radius R is 



/ 



ds-E = 2-KQ (1.96) 



where the 2it (instead of 47r) appears because the circumference of a circle is 
2-kR (and is analogous to AttR'^). Thus we find 

E{r) = ^ (1.97) 

and the 2D version of Poisson's equation is 

V • ^ = -VV = 27rQ 8'^{r). (1.99) 

Here ro is an arbitrary scale factor whose value is immaterial since it only shifts 
by a constant. 

We now see why the potential energy of interaction among a group of objects 
with charge m is 

UQ = m^^{-\u\z^- Zj\). (1.100) 

(Since z = {x + iy)/£ we are using ro = This explains the first term in 
eq. (1.91). 

To understand the second term notice that 

-V'\\z\' = -^=2npn (1.101) 
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where 

PB ^ (1.102) 

Eq. (1.101) can be interpreted as Poisson's equation and tells us that jlz]"^ 
represents the electrostatic potential of a constant charge density • Thus the 
second term in eq. (1.91) is the energy of charge m objects interacting with this 
negative background. 

Notice that 27r£^ is precisely the area containing one quantum of flux. Thus 
the background charge density is precisely B/^q, the density of flux in units of 
the flux quantum. 

The very long range forces in this fake plasma cost huge (fake) 'energy' unless 
the plasma is everywhere locally neutral (on length scales larger than the Debye 
screening length which in this case is comparable to the particle spacing). In 
order to be neutral, the density n of particles must obey 

nm + pB = (1.103) 

since each particle carries (fake) charge m. For our filled Landau level with 

m = 1, this is of course the correct answer for the density since every single- 
particle state is occupied and there is one state per quantum of flux. 

We again emphasize that the energy of the fake plasma has nothing to do 
with the quantum Hamiltonian and the true energy. The plasma analogy is 
merely a statement about this particular choice of wave function. It says that 
the square of the wave function is very small (because f/ciass is large) for config- 
urations in which the density deviates even a small amount from l/(27r£^). The 
electrons can in principle be found anywhere, but the overwhelming probability 
is that they are found in a configuration which is locally random (liquid-like) 
but with negligible density fluctuations on long length scales. We will discuss 
the nature of the typical configurations again further below in connection with 
fig. (1.12). 

When the fractional quantum Hall effect was discovered, Robert Laughlin 
realized that one could write down a many-body variational wave function at 
filling factor ly = 1/m by simply taking the mth power of the polynomial that 
describes the filled Landau level 

N 

fN[z] = l[('i-zjr- (1-105) 

In order for this to remain analytic, m must be an integer. To preserve the 
antisymmetry m must be restricted to the odd integers. In the plasma analogy 
the particles now have fake charge m (rather than unity) and the density of 
electrons is n = ^ so the Landau level filling factor v = ^ = i, i, i, 
etc. (Later on, other wave functions were developed to describe more general 
states in the hierarchy of rational fractional filling factors at which quantized 
Hall plateaus were observed [3,4,6,8,9].) 
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The Laughlin wave function naturally builds in good correlations among 
the electrons because each particle sees an m-fold zero at the positions of all 
the other particles. The wave function vanishes extremely rapidly if any two 
particles approach each other, and this helps minimize the expectation value of 
the Coulomb energy. 

Since the kinetic energy is fixed we need only concern ourselves with the 
expectation value of the potential energy for this variational wave function. 
Despite the fact that there are no adjustable variational parameters (other than 
m which controls the density) the Laughlin wave functions have proven to be 
very nearly exact for almost any realistic form of repulsive interaction. To 
understand how this can be so, it is instructive to consider a model for which 
this wave function actually is the exact ground state. Notice that the form of 
the wave function guarantees that every pair of particles has relative angular 
momentum greater than or equal to m. One should not make the mistake of 
thinking that every pair has relative angular momentum precisely equal to m. 
This would require the spatial separation between particles to be very nearly 
the same for every pair, which is of course impossible. 

Suppose that we write the Hamiltonian in terms of the Haldane pseudopo- 
tentials 

oo 

V=J2 E^-' Prn'iij) (1-106) 
m'=0 i<j 

where Pm{ij) is the projection operator which selects out states in which par- 
ticles i and j have relative angular momentum m. If Pm'iij) and P^"{jk) 
commuted with each other things would be simple to solve, but this is not 
the case. However if wc; (;onsi(k;r the; case of a 'hard-core potential' defined by 
Vm' = for to' > m, then clearly the TOth Laughlin state is an exact, zero energy 
eigenstate 

ViJm[z]=0. (1.107) 
This follows from the fact that 



Pm'{ij)^m = (L108) 

for any m' < m since every pair has relative angular momentum of at least to. 

Because the relative angular momentum of a pair can change only in discrete 
(even integer) units, it turns out that this hard core model has an excitation 
gap. For example for to = 3, any excitation out of the Laughlin ground state 
necessarily weakens the nearly ideal correlations by forcing at least one pair of 
particles to have relative angular momentum 1 instead of 3 (or larger). This 
costs an excitation energy of order vi . 

This excitation gap is essential to the existence of dissipationless {axx = 
Pxx = 0) current flow. In addition this gap means that the Laughlin state is 
stable against perturbations. Thus the difference between the Haldane pscu- 
dopotentials for the Coulomb interaction and the pseudopotentials for the 
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Figure 1.12: Comparison of typical configurations for a completely uncorrelated 
(Poissoii) distribution of 1000 particles (left panel) to the distribution given by 
the Laughlin wave function for m = 3 (right panel). The latter is a snapshot 
taken during a Monte Carlo simulation of the distribution. The Monte Carlo 
procedure consists of proposing a random trial move of one of the particles to 
a new position. If this move increases the value of j^Pp it is always accepted. If 
the move decreases the value of by a factor p, then the move is accepted 
with probability p. After equilibration of the plasma by a large number of such 
moves one finds that the configurations generated are distributed according to 
(After R. B. Laughlin, Chap. 7 in [3].) 

hard core model can be treated as a small perturbation (relative to the excita- 
tion gap). Numerical studies show that for realistic pscudopotcntials the overlap 
between the true ground state and the Laughlin state is extremely good. 

To get a better understanding of the correlations built into the Laughlin 
wave function it is useful to consider the snapshot in fig. (1.12) which shows a 
typical configuration of particles in the Laughlin ground state (obtained from 
a Monte Carlo sampling of |V'|^) compared to a random (Poisson) distribution. 
Focussing first on the large scale features wc sec that density fluctuations at 
long wavelengths are severely suppressed in the Laughlin state. This is easily 
understood in terms of the plasma analogy and the desire for local neutrality. A 
simple estimate for the density fluctuations p^- at wave vector q can be obtained 
by noting that the fake plasma potential energy can be written (ignoring a 
constant associated with self-interactions being included) 

f/class = ^ P^P-? (1-109) 

where is the area of the system and ^ is the Fourier transform of the 
logarithmic potential (easily derived from (— In (r)) = — 27r 5^{f) ). At long 
wavelengths (g^ <C n) it is legitimate to treat as a collective coordinate of an 
elastic continuum. The distribution e~'^^'='^'"' of these coordinates is a gaussian 
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Figure 1.13: Plot of the two-point correlation function h{r) = 1 — g{r) for the 
Laughlin plasma with = m = 3 (left panel) and m = 5 (right panel). Notice 

that, unlike the result for m = 1 given in cq. (1.112), g{r) exhibits the oscillatory 
behavior characteristic of a strongly coupled plasma with short-range solid-like 
local order. 



and so obeys (taking into account the fact that = [pq)*) 

We clearly see that the long-range (fake) forces in the (fake) plasma strongly 
suppress long wavelength density fluctuations. We will return more to this 
point later when we study collective density wave excitations above the Laughlin 
ground state. 

The density fluctuations on short length scales are best studied in real space. 
The radial correlation g{r) function is a convenient object to consider. g{r) tells 
us the density at r given that there is a particle at the origin 

9{r) = ^^^5^ / "^^^^ ■ • • / ^'-"^ ^' • • • , zn)\^ (1.111) 

where Z = (V'lV')) 'T' is the density (assumed uniform) and the remaining factors 
account for all the different pairs of particles that could contribute. The factors 
of density are included in the denominator so that lim^^oo 9{r) = 1- 

Because the m = 1 state is a single Slater determinant g{z) can be computed 
exactly 

5(^) = 1 - e-^l^l'. (1-112) 

Fig. (1.13) shows numerical estimates of h{r) = 1 — g{r) for the cases m = 3 
and 5. Notice that for the p = 1/m state g{z) ~ j^p™ for small distances. 
Because of the strong suppression of density fluctuations at long wavelengths, 
g{z) converges exponentially rapidly to unity at large distances. For m > 1, g 
develops oscillations indicative of solid-like correlations and, the plasma actually 
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freezes^ at m i=a 65. The Coulomb interaction energy can be expressed in terms 

of g{z) as^° 



^jd'z^^ [g[z) - 1] (1.113) 



where the ( — 1) term accounts for the neutrahzing background and e is the 
dielectric constant of the host semiconductor. We can interpret g(z) — 1 as the 
density of the 'exchange-correlation hole' surrounding each particle. 
The correlation energies per particle for m = 3 and 5 are [27] 



and 



1 mv\^^) 



1 mv\^B) 

N (V'slV's) 



-0.4100 ±0.0001 (1-114) 



= -0.3277 ±0.0002 (1.115) 



in units of e^/el which is « 161 K for e = 12.8 (the value in GaAs), B = lOT. 
For the filled Landau level (m = 1) the exchange energy is — \/f' as can be seen 
from eqs. (1.112) and (1.113). 



Exercise 1.15 Find the radial distribution function for a one- dimensional 

spinless free electron gas of density n by writing the ground state wave function 
as a single Slater determinant and then integrating out all but two of the 
coordinates. Use this first quantization method even if you already know how 
to do this calculation using second quantization. Hint: Take advantage of the 
following representation of the determinant of a N x N matrix M in terms 
of permutations P of N objects. 

N 



DetM = ^(-l)^n^i^r 



Exercise 1.16 Using the same method derive eq. (1.112). 



^That is, Monte Carlo simulation of shows that the particles are most likely to be found 
in a crystalline configuration which breaks translation symmetry. Again we emphasize that 
this is a statement about the Laughlin variational wave function, not necessarily a statement 
about what the electrons actually do. It turns out that for m > ~ 7 the Laughlin wave 
function is no longer the best variational wave function. One can write down wave functions 
describing Wigner crystal states which have lower variational energy than the Laughlin liquid. 

^"This expression assumes a strictly zero thickness electron gas. Otherwise one must replace 
-7-r by — f^°° ds— i^MJ= where F is the wavefunction factor describing the quantum well 
bound state. 
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1.7 Neutral Collective Excitations 

So far we have studied one particular variational wave function and found that 
it has good correlations built into it as graphically illustrated in Fig. 1.12. To 
further bolster the case that this wave function captures the physics of the 
fractional Hall effect we must now demonstrate that there is finite energy cost 
to produce excitations above this ground state. In this section we will study the 
neutral collective excitations. We will examine the charged excitations in the 
next section. 

It turns out that the neutral excitations are phonon-like excitations similar 
to those in solids and in superfluid helium. We can therefore use a simple 
modification of Feynman's theory of the excitations in superfluid helium [28,29]. 

By way of introduction let us start with the simple harmonic oscillator. The 
ground state is of the form 

M^) (1-116) 

Suppose we did not know the excited state and tried to make a variational 
ansatz for it. Normally we think of the variational method as applying only to 
ground states. However it is not hard to see that the flrst excited state energy 
is given by 

provided that we do the minimization over the set of states ip which are con- 
strained to be orthogonal to the ground state V-'o- One simple way to produce 
a variational state which is automatically orthogonal to the ground state is to 
change the parity by multiplying by the first power of the coordinate 

V-ilx) -xe""^'. (1.118) 

Variation with respect to a of course leads (in this special case) to the exact 
first excited state. 

With this background let us now consider the case of phonons in superfluid 
^He. Feynman argued that because of the Bose statistics of the particles, there 
are no low-lying single-particle excitations. This is in stark contrast to a fermi 
gas which has a high density of low-lying excitations around the fermi surface. 
Feynman argued that the only low-lying excitations in ^He are collective density 
oscillations that are well-described by the following family of variational wave 
functions (that has no adjustable parameters) labeled by the wave vector 

V'E = -^Pfc*o (1.119) 



where $o is the exact ground state and 

JV 
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Figure 1.14: (a) Configuration of particles in which the Fourier transform of 
the density at wave vector k is non-zero, (b) The Fourier amphtude will have a 
similar magnitude for this configuration but a different phase. 

is the Fourier transform of the density. The physical picture behind this is that 
at long wavelengths the fluid acts like an elastic continuum and can be treated 
as a generalized oscillator normal-mode coordinate. In this sense eq. (1.119) is 
then analogous to eq. (1.118). To see that tpj^ is orthogonal to the ground state 
we simply note that 

= -^Jd'Re-''-^{^o\p{r)\^o). (1.121) 

where 

N 

p{r)^Y.^\r,-R) (1.122) 
i=i 

is the density operator. If $o describes a translationally invariant liquid ground 
state then the Fourier transform of the mean density vanishes for k ^ 0. 

There are several reasons why ipj: is a good variational wave function, es- 
pecially for small k. First, it contains the ground state as a factor. Hence it 
contains all the special correlations built into the ground state to make sure that 
the particles avoid close approaches to each other without paying a high price in 
kinetic energy. Second, tpj: builds in the features we expect on physical grounds 
for a density wave. To see this, consider evaluating for a configuration of the 
particles like that shown in fig. (1.14a) which has a density modulation at wave 
vector k. This is not a configuration that maximizes |3>oP, but as long as the 
density modulation is not too large and the particles avoid close approaches, 
|$oP will not fall too far below its maximum value. More importantly, 
will be much larger than it would for a more nearly uniform distribution of posi- 
tions. As a result IV'gP will be large and this will be a likely configuration of the 
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particles in the excited state. For a configuration like that in fig. (1.14b), the 
phase of will shift but will have the same magnitude. This is analogous 
to the parity change in the harmonic oscillator example. Because all different 
phases of the density wave are equally likely, has a mean density which is 
uniform (translationally invariant). 

To proceed with the calculation of the variational estimate for the excitation 
energy A(A;) of the density wave state we write 

A(.) = M (1.123) 

where 

f{k)^(i;^\{H-EoM^), (1.124) 
with Eq being the exact ground state energy and 

s{k) ^ (^glV',-) = ^ ($o|pipg|$o). (1.125) 

We see that the norm of the variational state s{k) turns out to be the static 
structure factor of the ground state. It is a measure of the mean square density 
fluctuations at wave vector k. Continuing the harmonic oscillator analogy, we 
can view this as a measure of the zero-point fluctuations of the normal-mode 
oscillator coordinate pj:. For superfluid ''He s{k) can be directly measured by 
neutron scattering and can also be computed theoretically using quantum Monte 
Carlo methods [30]. We will return to this point shortly. 



Exercise 1.17 Show that for a uniform, liquid state of density n, the static 
structure factor is related to the Fourier transform of the radial distribution 
function by 



s{k) = N (5^ 3 + l + njd\ e"^-^ [g{r) - 1] 



The numerator in eq. (1.124) is called the oscillator strength and can be 
written 

/(fc) = ^ {Mp\[H,P^]\^^). (1.126) 

For uniform systems with parity symmetry we can write this as a double com- 
mutator 

/(fc) = ^ (*o \pUh,P^]\ (1.127) 



from which we can derive the justifiably famous oscillator strength sum rule 
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where M is the (band) mass of the particles. Remarkably (and conveniently) 
this is a universal result independent of the form of the interaction potential 
between the particles. This follows from the fact that only the kinetic energy 
part of the Hamiltonian fails to commute with the density. 



Exercise 1.18 Derive eg. (1.127) and then eq. (1.128) from eq. (1.126) for 

a system of interacting particles. 



We thus arrive at the Feynman-Bijl formula for the collective mode excitation 
energy 

We can interpret the first term as the energy cost if a single particle (initially at 
rest) were to absorb all the momentum and the second term is a renormalization 
factor describing momentum (and position) correlations among the particles. 
One of the remarkable features of the Feynman-Bijl formula is that it manages 
to express a dynamical quantity A(A:), which is a property of the excited state 
spectrum, solely in terms of a static property of the ground state, namely s{k). 
This is a very powerful and useful approximation. 

Returning to eq. (1.119) we see that ip^ describes a linear superposition of 

states in which one single particle has had its momentum boosted by hk. We 
do not know which one however. The summation in cq. (1.120) tells us that it 
is equally likely to be particle 1 or particle 2 or ... , etc. This state should not 
be confused with the state in which boost is applied to particle 1 and particle 
2 and ... , etc. This state is described by a product 





= I I I ^ ' I *o (1.130) 
which can be rewritten 

$,^ = exp|i7Vfc. l^^l^r.J j$o (1-131) 

showing that this is an exact energy eigenstatc (with energy N \^)m which 

the center of mass momentum has been boosted by NTik. 

In superfluid ^He the structure factor vanishes linearly at small wave vectors 

s{k) ~ ^fc (1.132) 

so that A(A;) is linear as expected for a sound mode 

A(^) ={ii-^k (1.133) 



Later on in Eq. (1.137) we will express the oscillator strength in terms of a frequency 
integral. Strictly spealcing if this is integrated up to very high frequencies including interband 
transitions, then M is replaced by the bare electron mass. 
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Figure 1.15: Schematic illustration of the phonon dispersion in superfluid liquid 
^He. For small wave vectors the dispersion is linear, as is expected for a gapless 
Goldstonc mode. The roton minimum due to the peak in the static structure 
factor occurs at a wave vector k of approximately 20 in units of inverse A. The 
roton energy is approximately 10 in units of Kelvins. 



from which we see that the sound velocity is given by 

'-mi 

This phonon mode should not be confused with the ordinary hydrodynamic 
sound mode in classical fluids. The latter occurs in a collision dominated regime 
uiT <i; 1 in which collision-induced pressure provides the restoring force. The 
phonon mode described here by is a low-lying eigenstate of the quantum 
Hamiltonian. 

At larger wave vectors there is a peak in the static structure factor caused 
by the solid- like oscillations in the radial distribution function g{r) similar to 
those shown in Fig. 1.13 for the Laughlin liquid. This peak in s{k) leads to the 
so-called roton minimum in A{k) as illustrated in fig. (1.15). 

To better understand the Feynman picture of the collective excited states 
recall that the dynamical structure factor is defined (at zero temperature) by 



t a/ H-Eo] 



$0^ • (1.135) 

The static structure factor is the zeroth frequency moment 

^^S{q,u;)= '^S{q,uj) (1.136) 

(with the second equality valid only at zero temperature). Similarly the oscil- 
lator strength in eq. (1.124) becomes (at zero temperature) 



f{q) 



j '^^^^1^'^)^! '^T^ S{q,uj). (1.137) 
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Thus we arrive at the result that the Feynman-Bijl formula can be rewritten 

w s huj S{q,uj) 

A(q) = :^a22_ ^^_^3g 

lo tSi<lM 

That is, A{q) is the mean excitation energy (weighted by the square of the 
density operator matrix element). Clearly the mean exceeds the minimum and 
so the estimate is variational as claimed. Feynman's approximation is equivalent 
to the assumption that only a single mode contributes any oscillator strength 
so that the zero-temperature dynamical structure factor contains only a single 
delta function peak 

S{q, u) = 2tt s{q) S (uj - ^ A(g)^ ■ (1.139) 

Notice that this approximate form satisfies both eq. (1.136) and eq. (1.137) 
provided that the collective mode energy A{q) obeys the Feynman-Bijl formula 
in eq. (1.129). 



Exercise 1.19 For a system with a homogeneous liquid ground state, the 
(linear response) static susceptibility of the density to a perturbation U = 
is defined by 

(p?) = x{q)v,-- (1-140) 

Using first order perturbation theory show that the static susceptibility is given 
in terms of the dynamical structure factor by 

Using the single mode approximation and the oscillator strength sum rule, 
derive an expression for the collective mode dispersion in terms ofx{q)- (Your 
answer should not involve the static structure factor. Note also that eq.(1.14.0) 
is not needed to produce the answer to this part. Just work with eq. (1.141)-) 



As we mentioned previously Feynman argued that in ^He the Bose sym- 
metry of the wave functions guarantees that unlike in Fermi systems, there is 
only a single low-lying mode, namely the phonon density mode. The paucity 
of low-energy single particle excitations in boson systems is what helps make 
them superfluid-there are no dissipative channels for the current to decay into. 
Despite the fact that the quantum Hall system is made up of fermions, the 
behavior is also reminiscent of superfluidity since the current flow is dissipa- 
tionless. Indeed, within the 'composite boson' picture, one views the FQHE 
ground state as a bose condensate [1,9,10]. Let us therefore blindly make the 
single-mode approximation and see what happens. 

From eq. (1.110) we see that the static structure factor for the mth Laughlin 
state is (for small wave vectors only) 
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where we have used L^/N = 2'jTi'^m. The Feynman-Bijl formula then yields^^ 

This predicts that there is an excitation gap that is independent of wave vector 
(for small q) and equal to the cyclotron energy. It is in fact correct that at 
long wavelengths the oscillator strength is dominated by transitions in which a 
single particle is excited from the n = to the n = 1 Landau level. Furthermore, 
Kohn's theorem guarantees that the mode energy is precisely huic- Eq. (1.143) 
was derived specifically for the Laughlin state, but it is actually quite general, 
applying to any translationally invariant liquid ground state. 

One might expect that the single mode approximation (SMA) will not work 
well in an ordinary Fermi gas due to the high density of excitations around the 
Fermi surface. Here however the Fermi surface has been destroyed by the 
magnetic field and the continuum of excitations with different kinetic energies 
has been turned into a set of discrete inter-Landau-level excitations, the lowest 
of which dominates the oscillator strength. 

For filling factor u = 1 the Pauli principle prevents any intra-level excitations 
and the excitation gap is in fact KlUc as predicted by the SMA. However for i/ < I 
there should exist intra-Landau-level excitations whose energy scale is set by the 
interaction scale e^/e£ rather than the kinetic energy scale Uuic- Indeed we can 
formally think of taking the band mass to zero {M — > 0) which would send 
hiVc oo while keeping e^/ei fixed. Unfortunately the SMA as it stands now 
is not very useful in this limit. What we need is a variational wave function 
that represents a density wave but is restricted to lie in the Hilbert space of the 
lowest Landau level. This can be formally accomplished by replacing eq. (1.119) 

by 

= Pk ^rn (1.144) 

where the overbar indicates that the density operator has been projected into 
the lowest Landau level. The details of how this is accomplished are presented 
in appendix A. 

The analog of eq. (1.123) is 

A(.) = ffl (1.145) 

where / and s are the projected oscillator strength and structure factor, respec- 
tively. As shown in appendix A 



s{k) = ^ (V™|4 PfelV'™) = s{k) - [l - e-il*=l 



= s{k)- s,=x{k). (1.146) 



^^We will continue to use the symbol M here for the band mass of the electrons to avoid 
confusion with the inverse filling factor m. 

^^This expectation is only partly correct however as one discovers when studying collective 
plasma oscillations in systems with long-range Coulomb forces. 
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This vanishes for the filled Landau level because the Pauli principle forbids all 
intra-Landau-level excitations. For the mth Laughlin state eq. (1.142) shows us 
that the leading term in s{k) for small k is ^k'^i'^. Putting this into eq. (1.146) 
we see that the leading behavior for s{k) is therefore quartic 

s{k) a{k£)'^ + . . . . (1.147) 

We can not compute the coefficient a without finding the fc* correction to 
eq. (1.142). It turns out that there exists a compressibility sum rule for the 
fake plasma from which we can obtain the exact result [29] 

a=^. (1.148) 

The projected oscillator strength is given by eq. (1.127) with the density 

operators replaced by their projections. In the case of ^He only the kinetic 
energy part of the Hamiltonian failed to commute with the density. It was for 
this reason that the oscillator strength came out to be a universal number related 
to the mass of the particles. Within the lowest Landau level however the kinetic 
energy is an irrelevant constant. Instead, after projection the density operators 
no longer commute with each other (see appendix A). It follows from these 
commutation relations that the projected oscillator strength is proportional to 
the strength of the interaction term. The leading small k behavior is [29] 

/(fc) = 6^(M)4 + ... (1.149) 

where 6 is a dimensionless constant that depends on the details of the interaction 
potential. The intra-Landau level excitation energy therefore has a finite gap 
at small k 

This is quite different from the case of superfluid ^He in which the mode is 
gapless. However like the case of the superfluid, this 'magnetophonon' mode 
has a 'magnetoroton' minimum at finite k as illustrated in fig. (1.16). The 
figure also shows results from numerical exact diagonalization studies which 
demonstrate that the single mode approximation is extremely accurate. Note 
that the magnetoroton minimum occurs close to the position of the smallest 
reciprocal lattice vector in the Wigner crystal of the same density. In the crystal 
the phonon frequency would go exactly to zero at this point. (Recall that in 
a crystal the phonon dispersion curves have the periodicity of the reciprocal 
lattice.) 

Because the oscillator strength is almost entirely in the cyclotron mode, the 
dipole matrix element for coupling the collective excitations to light is very 
small. They have however been observed in Raman scattering [33] and found 
to have an energy gap in excellent quantitative agreement with the single mode 
approximation. 
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Figure 1.16: Comparison of the single mode approximation (SMA) prediction of 
the collective mode energy for filling factors v = 1/3, 1/5, 1/7 (solid lines) with 
small-system numerical results for N particles. Crosses indicate the N = 1 ,v = 
1/3 spherical system, triangles indicate the iV = 6, i^=l/3 hexagonal miit cell 
system results of Haldane and Rezayi [31]. Solid dots are for = 9, = 1/3 and 
N = 7,v = 1/5 spherical system calculations of Fano et al. [32] Arrows at the 
top indicate the magnitude of the reciprocal lattice vector of the Wigner crystal 
at the corresponding filling factor. Notice that unlike the phonon collective 
mode in superfluid helium shown in fig. (1.15), the mode here is gapped. 
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Finally wc remark that these collective excitations are characterized by a 
well-defined wave vector k despite the presence of the strong magnetic field. 
This is only possible because they are charge neutral which allows one to define 
a gauge invariant conserved momentum [34]. 



1.8 Charged Excitations 

Except for the fact that they are gapped, the neutral magnetophonon excitations 
are closely analogous to the phonon excitations in superfluid ^He. We further 
pursue this analogy with a search for the analog of vortices in superfluid films. 
A vortex is a topological defect which is the quantum version of the familiar 
whirlpool. A reasonably good variational wave function for a vortex in a two- 
dimensional film of *He is 

V'^ = I n / (l^J- - ^l) e^''^^^^-^^ i $0. (1.151) 



Here 6 is the azimuthal angle that the particle's position makes relative to R, 
the location of the vortex center. The function / vanishes as f approaches R 
and goes to unity far away. The choice of sign in the phase determines whether 
the vortex is right or left handed. 

The interpretation of this wave function is the following. The vortex is 
a topological defect because if any particle is dragged around a closed loop 
surrounding R, the phase of the wave function winds by ±27r. This phase 
gradient means that current is circulating around the core. Consider a large 
circle of radius ^ centered on R. The phase change of 27r around the circle 
occurs in a distance 27r^ so the local gradient seen by every particle is 9/^. 
Recalling eq. (1.131) we see that locally the center of mass momentum has 
been boosted by ±| ^ so that the current density of the whirlpool falls off 
inversely with distance from the core.^^ Near the core / falls to zero because 
of the 'centrifugal barrier' associated with this circulation. In a more accurate 
variational wave function the core would be treated slightly differently but the 
asymptotic large distance behavior would be unchanged. 

What is the analog of all this for the lowest Landau level? For ^+ we see 
that every particle has its angular momentum boosted by one unit. In the lowest 
Landau level analyticity (in the symmetric gauge) requires us to replace e'^ by 
z = X + iy. Thus we are led to the Laughlin 'quasi-hole' wave function 

N 

^tM = Y[{zj-Z)^m[z] (1.152) 



^*This slow algebraic decay of the current density means that the total kinetic energy of a 
single vortex diverges logarithmically with the size of the system. This in turn leads to the 
Kosterlitz Thouless phase transition in which pairs of vortices bind together below a critical 
temperature. As we will see below there is no corresponding finite temperature transition in 
a quantum Hall system. 
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where Z is a complex miniber denoting the position of the vortex and ■i/'m is the 
Laughhn wave function at fiUing factor v = 1/m. The corresponding antivortex 
('quasi-electron' state) involves Zj suitably projected (as discussed in App. A.): 



N 



V-zN = n ( 2^ - ^* ) '^'"W (1-153) 



3 = 1 



where as usual the derivatives act only on the polynomial part of tprn- All 
these derivatives make ■(/'" somewhat difficult to work with. We will therefore 
concentrate on the quasi- hole state V'"'"- The origin of the names quasi- hole and 
quasi-electron will become clear shortly. 

Unlike the case of a superfluid film, the presence of the vector potential 
allows these vortices to cost only a finite energy to produce and hence the 
electrical dissipation is always finite at any non-zero temperature. There is no 
finite temperature transition into a superfluid state as in the Kosterlitz Thouless 
transition. From a field theoretic point of view, this is closely analogous to the 
Higg's mechanism [1]. 

Just as in our study of the Laughlin wave function, it is very useful to see 
how the plasma analogy works for the quasi-hole state 

|V;+|2 = e-'^^'=— e-/^^ (1.154) 

where J/dass is given by eq. (1.91), (3 = 2/m as before and 

N 

l/^m^(-ln|z,-Z|). (1.155) 

Thus we have the classical statistical mechanics of a one-component plasma of 
(fake) charge m, objcc;ts seeing a neutralizing jellium background phis a new po- 
tential energy V representing the interaction of these objects with an 'impurity' 
located at Z and having unit charge. 

Recall that the chief desire of the plasma is to maintain charge neutrality. 
Hence the plasma particles will be repelled from Z. Because the plasma particles 
have fake charge m, the screening cloud will have to have a net reduction of 
1/m particles to screen the impurity. But this means that the quasi-hole has 
fractional fermion number! The (true) physical charge of the object is a fraction 
of the elementary charge 

q* = -. (1.156) 
m 

This is very strange! How can we possibly have an elementary excitation 
carrying fractional charge in a system made up entirely of electrons? To under- 
stand this let us consider an example of another quantum system that seems to 
have fractional charge, but in reality doesn't. Imagine three protons arranged 
in an equilateral triangle as shown in fig. (1.17). Let there be one electron in 
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IS 

Figure 1.17: Illustration of an electron tunneling among the IS orbitals of three 
protons. The tunneling is exponentially slow for large separations which leads 
to only exponentially small lifting of what would otherwise be a three-fold de- 
generate ground state. 



the system. In the spirit of the tight-binding model wc consider only the IS 
orbital on each of the three 'lattice sites'. The Bloch states are 

1 ^ 

i^k = -^J2^'''\j) (1.157) 

where \j) is the IS orbital for the jth atom. The equilateral triangle is like 
a linear system of length 3 with periodic boundary conditions. Hence the al- 
lowed values of the wavevector are {ka = ^a; a = — l,0,-|-l}. The energy 
eigenvalues are 

efe^ = -Eis - 2 J cosfca (1.158) 

where E'ls is the isolated atom energy and — J is the hopping matrix element 
related to the orbital overlap and is exponentially small for large separations of 
the atoms. 

The projection operator that measures whether or not the particle is on site 
n is 

P„ = \n){n\. (1.159) 
Its expectation value in any of the three eigenstates is 

{^kJPn\i^kJ = l. (1.160) 

This equation simply reflects the fact that as the particle tunnels from site to 
site it is equally likely to be found on any site. Hence it will, on average, be 
found on a particular site n only 1/3 of the time. The average electron number 
per site is thus 1/3. This however is a trivial example because the value of 
the measured charge is always an integer. Two-thirds of the time we measure 
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zero and one third of the time we measure unity. This means that the charge 
fluctuates. One measure of the fluctuations is 



which shows that the fluctuations are larger than the mean value. This result 
is most easily obtained by noting P^^ — Pn- 

A characteristic feature of this 'imposter' fractional charge ^ that guaran- 
tees that it fluctuates is the existence in the spectrum of the Hamiltonian of 
a set of TO nearly degenerate states. (In our toy example here, to ~ 3.) The 
characteristic time scale for the charge fluctuations is t ~ h/Ae where Ae is the 
energy splitting of the quasi-degenerate manifold of states. In our tight-binding 
example t ^ h/J is the characteristic time it takes an electron to tunnel from 
the IS orbital on one site to the next. As the separation between the sites in- 
creases this tunneling time grows exponentially large and the charge fluctuations 
become exponentially slow and thus easy to detect. 

In a certain precise sense, the fractional charge of the Laughlin quasiparticles 
behaves very diff'erently from this. An electron added at low energies to a = 
1/3 quantum Hall fluid breaks up into three charge 1/3 Laughlin quasiparticles. 
These quasiparticles can move arbitrarily far apart from each other^^ and yet 
no quasi-degenerate manifold of states appears. The excitation gap to the first 
excited state remains finite. The only degeneracy is that associated with the 
positions of the quasiparticles. If we imagine that there are three impurity 
potentials that pin down the positions of the three quasiparticles, then the state 
of the system is uniquely specified. Becaiisc there is no qiiasidegcneracy, we 
do not have to specify any more information other than the positions of the 
quasiparticles. Hence in a deep sense, they are true elementary particles whose 
fractional charge is a sharp quantum observable. 

Of course, since the system is made up only of electrons, if we capture the 
charges in some region in a box, we will always get an integer number of electrons 
inside the box. However in order to close the box we have to locally destroy the 
Laughlin state. This will cost (at a minimum) the excitation gap. This may not 
seem important since the gap is small — only a few Kelvin or so. But imagine 
that the gap were an MeV or a GeV. Then we would have to build a particle 
accelerator to 'close the box' and probe the fluctuations in the charge. These 
fluctuations would be analogous to the ones seen in quantum electrodynamics 
at energies above 2raeC? where electron-positron pairs are produced during the 
measurement of charge form factors by means of a scattering experiment. 

Put another way, the charge of the Laughlin quasiparticle fluctuates but 
only at high frequencies ^ /S./h. If this frequency (which is ^ 50GHz) is higher 
than the frequency response limit of our voltage probes, we will see no charge 
fluctuations. We can formalize this by writing a modified projection operator 
[35] for the charge on some site n by 




(1.161) 




(1.162) 



Recall that unlike the case of vortices in superfluids, these objects are unconfined. 
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where P„ = |n) {n\ as before and 

P^^'> = 9{n - H + Eo) (1.163) 

is the operator that projects onto the subset of eigenstates with excitation en- 
ergies less than n. P^'^ thus represents a mcasiircmcnt with a high-frequency 
cutoff built in to represent the finite bandwidth of the detector. Returning to 
our tight-binding example, consider the situation where J is large enough that 
the excitation gap A = (l — cos ^) J exceeds the cutoff f2. Then 

+1 

a=-l 

= IV'feo)(V'fcol (1-164) 
is simply a projector on the ground state. In this case 

^f) = IV'feo)^(^feol (1-165) 



and 

>feo 



[^f ¥ V'fco) - (V'fcol^f Vfco)' = 0. (1.166) 



Tlie charge fluctuations in the ground state are then zero (as measured by the 
finite bandwidth detector). 

The argument for the Laughlin quasiparticles is similar. We again emphasize 
that one can not think of a single charge tunneling among three sites because the 
excitation gap remains finite no matter how far apart the quasiparticle sites are 
located. This is possible only because it is a correlated many-particle system. 

To gain a better understanding of fractional charge it is useful to compare 
this situation to that in high energy physics. In that field of study one knows 
the physics at low energies — this is just the phenomena of our everyday world. 
The goal is to study the high energy (short length scale) limit to see where this 
low energy physics comes from. What force laws lead to our world? Probing 
the proton with high energy electrons we can temporarily break it up into three 
fractionally charged quarks, for example. 

Condensed matter physics in a sense does the reverse. We know the phe- 
nomena at 'high' energies (i.e. room temperature) and we would like to see how 
the known dynamics (Coulomb's law and non-relativistic quantum mechanics) 
leads to unknown and surprising collective effects at low temperatures and long 
length scales. The analog of the particle accelerator is the dilution refrigerator. 

To further understand Laughlin quasiparticles consider the point of view 
of 'flatland' physicists living in the cold, two-dimensional world of a = 1/3 
quantum Hall sample. As far as the flatlanders are concerned the 'vacuum' 
(the Laughlin liquid) is completely inert and featureless. They discover however 
that the universe is not completely empty. There are a few elementary particles 
around, all having the same charge q. The flatland equivalent of Benjamin 
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Franklin chooses a unit of charge which not only makes q negative but gives it 
the fractional value —1/3. For some reason the Flatlanders go along with this. 

Flatland cosmologists theorize that these objects are 'cosmic strings', topo- 
logical defects left over from the 'big cool down' that followed the creation of 
the universe. Flatland experimentalists call for the creation of a national ac- 
celerator facility which will reach the unprecedented energy scale of 10 Kelvin. 
With great effort and expense this energy scale is reached and the accelerator 
is used to smash together three charged particles. To the astonishment of the 
entire world a new short-lived particle is temporarily created with the bizarre 
property of having integer charge! 

There is another way to see that the Laughlin quasiparticles carry fractional 
charge which is useful to understand because it shows the deep connection be- 
tween the sharp fractional charge and the sharp quantization of the Hall con- 
ductivity. Imagine piercing the sample with an infinitely thin magnetic solenoid 
as shown in fig. (1.18) and slowly increasing the magnetic flux $ from to 
$Q = ^ the quantum of flux. Because of the existence of a finite excitation gap 
A the process is adiabatic and reversible if performed slowly on a time scale 
long compared to h/A. 

Faraday's law tells us that the changing flux induces an electric field obeying 

.dr-E^ — ^ (1.167) 

where V is any contour surrounding the flux tube. Because the electric field 
contains only Fourier components at frequencies lo obeying ?ia> < A, there is no 
dissipation and Gxx = fyj/ = Pxx = Pyy = 0. The electric field induces a current 
density obeying 

E = PxyJxz (1.168) 
so that 

Pxy (f J-{zxdf) = --^. (1.169) 
Jr c dt 

The integral on the LHS represents the total current flowing into the region 
enclosed by the contour. Thus the charge inside this region obeys 

After one quantum of flux has been added the flnal charge is 

<5 = CTxy^O = CTxy. (1-171) 

2 

Thus on the quantized Hall plateau at filling factor i/ where axy = ^ we have 
the result 

Q = ve. (1.172) 
Reversing the sign of the added fiux would reverse the sign of the charge. 



The Quantum Hall Effect 



53 



0(t) 



AAA 




E(t) \ 




Figure 1.18: Construction of a Laughlin quasiparticlc by adiabatically threading 
flux through a point in the sample. Faraday induction gives an azimuthal 
electric field E{t) which in turn produces a radial current J{t). For each quan- 
tum of flux added, charge ve flows into (or out of) the region due to the quan- 
tized Hall conductivity ue'^/h. A flux tube containing an integer number of flux 
quanta is invisible to the particles (since the Aharanov phase shift is an integer 
multiple of 27r) and so can be removed by a singular gauge transformation. 
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The final step in the argument is to note that an infinitesimal tube containing 
a quantum of flux is invisible to the particles. This is because the Aharonov- 
Bohm phase factor for traveling around the flux tube is unity. 

exp (f SA- df] = 6^2^^* = 1. (1.173) 



he 



r 



Here 6A is the additional vector potential due to the solenoid. Assuming the 
flux tube is located at the origin and making the gauge choice 

5A = <i>o^, (1.174) 

one can see by direct substitution into the Schrodinger equation that the only 
effect of the quantized flux tube is to change the phase of the wave function by 

The removal of a quantized flux tube is thus a 'singular gauge cliange' which 
has no physical effect. 

Let us reiterate. Adiabatic insertion of a flux quantum changes the state of 
the system by pulling in (or pushing out) a (fractionally) quantized amount of 
charge. Once the flux tube contains a quantum of flux it effectively becomes 
invisible to the electrons and can be removed by means of a singular gauge 
transformation. 

Because the excitation gap is preserved during the adiabatic addition of the 
flux, the state of the system is fully specifled by the position of the resulting 
quasiparticle. As discussed before there are no low-lying quasi-degenerate states. 
This version of the argument highlights the essential importance of the fact that 
<^xx = and axy is quantized. The existence of the fractionally quantized Hall 
transport coefficients guarantees the existence of fractionally charged elementary 
excitations 

These fractionally charged objects have been observed directly by using an 
ultrasensitive electrometer made from a quantum dot [36] and by the reduced 
shot noise which they produce when they carry current [37]. 

Because the Laughlin quasiparticles are discrete objects they cost a non-zero 
(but finite) energy to produce. Since they are charged they can be thermally 
excited only in neutral pairs. The charge excitation gap is therefore 

Ae = A+ + A_ (1.176) 

where A± is the vortex/antivortex (quasielectron/quasihole) excitation energy. 

In the presence of a transport current these thermally excited charges can move 
under the influence of the Hall electric fleld and dissipate energy. The resulting 
resistivity has the Arrhenius form 



e 



7^ e-'^^^/^ (1.177) 
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where 7 is a dimcnsionlcss constant of order unity. Note that the law of mass 
action tells us that the activation energy is Ac/2 not Ag since the charges are 
excited in pairs. There is a close analogy between the dissipation described here 
and the flux flow resistance caused by vortices in a superconducting film. 

Theoretical estimates of A^ are in good agreement with experimental values 
determined from transport measurements [38]. Typical values of Ac are only a 
few percent of e? j el and hence no larger than a few Kelvin. In a superfiuid time- 
reversal symmetry guarantees that vortices and antivortices have equal energies. 
The lack of time reversal symmetry here means that A_|_ and A_ can be quite 
different. Consider for example the hard-core model for which the Laughlin 
wave function is an exact zero energy ground state as shown in eq. (1.107). 
Equation (1.152) shows that the quasihole state contains Vm as a factor and 
hence is also an exact zero energy eigenstate for the hard-core interaction. Thus 
the quasihole costs zero energy. On the other hand eq. (1.153) tells us that 
the derivatives reduce the degree of homogeneity of the Laughlin polynomial 
and therefore the energy of the quasiclectron must be non-zero in the hard-core 
model. At filling factor v = 1/m this asymmetry has no particular significance 
since the quasiparticles must be excited in pairs. 

Consider now what happens when the magnetic field is increased slightly 
or the particle number is decreased slightly so that the filling factor is slightly 
smaller than 1/m. The lowest energy way to accommodate this is to inject 
m quasiholes into the Laughlin state for each electron that is removed (or for 
each m^Q of fiux that is added). The system energy (ignoring disorder and 
interactions in the dilute gas of quasiparticles) is 

E+=Em- 5N mA+ (1.178) 

where Em is the Laughlin ground state energy and —SN is the number of added 
holes. Conversely for filling factors slightly greater than 1 /m the energy is (with 
+5N being the number of added electrons) 

E_=Em + SNmA_. (1.179) 

This is illustrated in fig. (1.19). The slope of the lines in the figure determines 
the chemical potential 

M± = 1^ = TmA±. (1.180) 

The chemical potential suffers a jump discontinuity of m(A+ -I- A_) = mAc 
just at filling factor /x = 1/m. This jump in the chemical potential is the 
signature of the charge excitation gap just as it is in a semiconductor or insulator. 
Notice that this form of the energy is very reminiscent of the energy of a type- 
II superconductor as a function of the applied magnetic field (which induces 
vortices and therefore has an energy cost AE ~ \B\). 

Recall that in order to have a quantized Hall plateau of finite width it is 
necessary to have disorder present. For the integer case we found that disorder 
localizes the excess electrons allowing the transport coefficients to not change 
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Figure 1.19: Energy cost for inserting 6N electrons into the Laughlin state near 
filling factor u = 1/m. The slope of the line is the chemical potential. Its 
discontinuity at v = 1/m measures the charge excitation gap. 



with the filling factor. Here it is the fractionally-charged quasiparticles that arc 
localized by the disorder. Just as in the integer case the disorder may fill in 
the gap in the density of states but the DC value of axx can remain zero because 
of the localization. Thus the fractional plateaus can have finite width. 

If the density of quasiparticles becomes too high they may delocalize and 
condense into a correlated Laughlin state of their own. This gives rise to a 
hierarchical family of Hall plateaus at rational fractional filling factors v = p/ q 
(generically with q odd due to the Pauli principle). There are several diflFerent 
but entirely equivalent ways of constructing and viewing this hierarchy which 
we will not delve into here [3,4,6]. 



1.9 FQHE Edge States 

We learned in our study of the integer QHE that gapless edge excitations exist 

even when the bulk has a large excitation gap. Because the bulk is incom- 
pressible the only gapless neutral excitations must be area-preserving shape 
distortions such as those illustrated for a disk geometry in fig. (1.20a). Because 
of the confining potential at the edges these shape distortions have a character- 
istic velocity produced by the E x B drift. It is possible to show that this view 
of the gapless neutral excitations is precisely equivalent to the usual Fermi gas 
particlc-holc pair excitations that we considered previously in our discussion of 
edge states. Recall that we argued that the contour line of the electrostatic 
potential separating the occupied from the empty states could be viewed as a 
real-space analog of the Fermi siirfacc (since position and momentum are equiv- 
alent in the Landau gauge). The charged excitations at the edge are simply 
ordinary electrons added or removed from the vicinity of the edge. 



^^Note again the essential importance of the fact that the objects are 'elementary particles'. 
That is, there are no residual degeneracies once the positions are pinned down. 
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Figure 1.20: Arca-prcscrving shape distortions of tlie ineompressible quantum 
Hall state, (a) IQHE Lauglilin liquid 'droplet' at v = 1. (b) FQHE annulus 
at u = 1/m formc;d by injeeting a large number n of flux quanta at the origin 
to create n quasiholes. There are thus two edge modes of opposite chirality. 
Changing n by one unit transfers fractional charge i>e from one edge to the other 
by expanding or shrinking the size of the central hole. Thus the edge modes 
have topological sectors labeled by the 'winding number' n and one can view the 
gapless edge excitations as a gas of fractionally charged Laughlin quasiparticles. 
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In the case of a fractional QHE state at = 1/m the bulk gap is caused 
by Coulomb correlations and is smaller but still finite. Again the only gapless 
excitations are area-preserving shape distortions. Now however the charge of 
each edge can be varied in units of e/m. Consider the annulus of Hall fluid 
shown in fig. (1.20b). The extension of the Laughlin wave function to this 
situation is 



This simply places a large number n :» 1 of quasiholes at the origin. Following 
the plasma analogy we see that this looks like a highly charged impurity at the 
origin which repels the plasma, producing the annulus shown in fig. (1.20b). 
Each time we increase n by one unit, the annulus expands. We can view this 
expansion as increasing the electron number at the outer edge by 1/m and 
reducing it by 1/m at the inner edge. (Thereby keeping the total electron 
number integral as it must be.) 

It is appropriate to view the Laughlin quasiparticles, which are gapped in 
the bulk, as being liberated at the edge. The gapless shape distortions in the 
Hall liquid are thus excitations in a 'gas' of fractionally charged quasiparticles. 
This fact produces a profound alteration in the tunneling density of states to 
inject an electron into the system. An electron which is suddenly added to 
an edge (by tunneling through a barrier from an external electrode) will have 
very high energy unless it breaks up into m Laughlin quasiparticles. This leads 
to an 'orthogonality catastrophe' which simply means that the probability for 
this process is smaller and smaller for final states of lower and lower energy. 
As a result the current-voltage characteristic for the tunnel junction becomes 
non-linear [17,39,40] 



For the filled Landau level m = 1 the quasiparticles have charge q = em = e 
and are ordinary electrons. Hence there is no orthogonality catastrophe and the 
I-V characteristic is linear as expected for an ordinary metallic tunnel junction. 
The non-linear tunneling for the m = 3 state is shown in fig. (1.21). 

1.10 Quantum Hall Ferromagnets 
1.10.1 Introduction 

Naively one might imagine that electrons in the QHE have their spin dynamics 
frozen out by the Zeeman splitting gjJL^B. In free space with g = 2 (neglecting 
QED corrections) the Zeeman splitting is exactly equal to the cyclotron splitting 
huc ~ 100 K as illustrated in fig. (1.22 a). Thus at low temperatures we would 
expect for filling factors < 1 all the spins would be fully aligned. It turns out 
however that this naive expectation is incorrect in GaAs for two reasons. First, 




(1.181) 



(1.182) 
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Figure 1.21: Non- linear current voltage response for tunneling an electron into 
a FQHE edge state. Because the electron must break up into m fractionally 
charged quasiparticles, there is an orthogonality catastrophe leading to a power- 
law density of states. The flattening at low currents is due to the finite tem- 
perature. The upper panel shows the f = 1/3 Hall plateau. The theory [17,39] 
works extremely well on the 1/3 quantized Hall plateau, but the unexpectedly 
smooth variation of the exponent with magnetic field away from the plateau 
shown in the lower panel is not yet fully understood. (After M. Grayson et al, 
Ref. [41]. 
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Figure 1.22: (a) Landau energy levels for an electron in free space. Numbers 
label the Landau levels and +(— ) refers to spin up (down). Since the g factor 
is 2, the Zeeman splitting is exactly equal to the Landau level spacing, Ticjc and 
there are extra degeneracies as indicated, (b) Same for an electron in GaAs. 
Because the effective mass is small and g « —0.4, the degeneracy is strongly 
lifted and the spin assignments are reversed. 
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the small effective mass (m* = 0.068) in the conduction band of GaAs increases 
the cyclotron energy by a factor of m/m* ~ 14. Second, spin-orbit scattering 
tumbles the spins around in a way which reduces their effective coupling to the 
external magnetic field by a factor of —5 making the g factor —0.4. The Zccman 
energy is thus some 70 times smaller than the cyclotron energy and typically 
has a value of about 2K, as indicated in fig. (1.22 b). 

This decoupling of the scales of the orbital and spin energies means that 
it is possible to be in a regime in which the orbital motion is fully quantized 
{ksT <^ fuvc) but the low-energy spin fluctuations are not completely frozen 
out (fceT ^ g*iiBB). The spin dynamics in this regime are extremely unusual 
and interesting because the system is an itinerant magnet with a quantized Hall 
coefficient. As we shall see, this leads to quite novel physical effects. 

The introduction of the spin degree of freedom means that we are dealing 
with the QHE in multicomponent systems. This subject has a long history going 
back to an early paper by Halperin [42] and has been reviewed extensively [4,43, 
44] . In addition to the spin degree of freedom there has been considerable recent 
interest in other multicomponent systems in which spin is replaced by a pseudo- 
spin representing the layer index in double well QHE systems or the electric 
subband index in wide single well systems. Experiments on these systems are 
discussed by Shayegan in this volume [45] and have also been reviewed in [44]. 

Our discussion will focus primarily on ferromagnetism near filling factor 
V = \. In the subsequent section we will address analogous effects for pseudo- 
spin degrees of freedom in multilayer systems. 

1.10.2 Coulomb Exchange 

We tend to think of the integer QHE as being associated with the gap due to 
the kinetic energy and ascribe importance to the Coulomb interaction only in 
the fractional QHE. However study of ferromagnetism near integer filling factor 
u = 1 has taught us that Coulomb interactions play an important role there as 
well [46]. 

Magnetism occurs not because of direct magnetic forces, but rather because 
of a combination of electrostatic forces and the Pauli principle. In a fully ferro- 
magnetically aligned state all the spins are parallel and hence the spin part of 
the wave function is exchange symmetric 



The spatial part $ of the wave function must therefore be fully antisymmetric 
and vanish when any two particles approac;li c;ac;h other. This means that each 
particle is surrounded by an 'exchange hole' which thus lowers the Coulomb 
energy per particle as shown in eq. (1.113). For filling factor u = 1 



This energy scale is two orders of magnitude larger than the Zeeman splitting 
and hence strongly stabilizes the ferromagnetic state. Indeed at v = 1 the 



|^) = $(zi,...,z^)|tTTTT...T). 



(1.183) 




(1.184) 
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ground state is spontaneously fully polarized at zero temperature even in the 
absence of the Zeeman term. Ordinary ferromagnets like iron are generally 
only partially polarized because of the extra kinetic energy cost of raising the 
fermi level for the majority carriers. Here however the kinetic energy has been 
quenched by the magnetic field and all states in the lowest Landau level are 
degenerate. For v = 1 the large gap to the next Landau level means that we 
know the spatial wave function <I> essentially exactly. It is simply the single 
Slater determinant representing the fully filled Landau level. That is, it is 
m = 1 Laughlin wave function. This simple circumstance makes this perhaps 
the world's best understood ferromagnet. 



1.10.3 Spin Wave Excitations 

It turns out that the low-lying 'magnon' (spin wave) excited states can also be 
obtained exactly. Before doing this for the QHE system let us remind ourselves 
how the calculation goes in the lattice Heisenberg model for A'' local moments 
in an insulating ferromagnet 

(ij) i 

= -J^|5f5; + i(5+57+5T5+)|-A^5| (1.185) 

The ground state for J > is the fully ferromagnetic state with total spin 

S = N/2. Let us choose our coordinates in spin space so that Sz = N/2. 
Because the spins are fully aligned the spin-flip terms in H are ineffective and 
(ignoring the Zeeman term) 

|TTT...T) = -^A^6|TTT...T) (1.186) 

where Ni, is the number of near- neighbor bonds and we have set h = 1. There 
are of course 25 -|- 1 = -|- 1 other states of the same total spin which will 
be degenerate in the absence of the Zeeman coupling. These are generated by 
successive applications of the total spin lowering operator 

N 

s- ^ J2^7 (1-187) 
j=l 

5-ITTT.-.T) = UTT.-.T) + ITiT..-T) 

+ I TTi...T) + ... (1.188) 

It is not hard to show that the one-magnon excited states are created by a 
closely related operator 

JV 

SZ ^ J2 e"'''^' Sr (1.189) 
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where g lies inside the Brillouin zone and is the magnon wave vector. Denote 

these states by 

m - SZ iV^o) (1.190) 

where lipa) is the ground state. Because there is one flipped spin in these states 
the transverse part of the Hcisenberg interaction is able to move the flipped spin 
from one site to a neighboring site 

J ^ 

-2 5IE^""-''^^^.^IV'o> (1.191) 

s J=i 

H\i;^) = {Eo + e^)\i;^) (1.192) 

where z is the coordination number, S is summed over near neighbor lattice 
vectors and the magnon energy is 

e,^^ |l-^E^""-'|+^ (1-193) 
For small q the dispersion is quadratic and for a 2D square lattice 

e,~-^q' + A (1.194) 

where a is the lattice constant. 

This is very different from the result for the antiferromagnet which has a 
linearly dispersing collective mode. There the ground and excited states can 
only be approximately determined because the ground state does not have all 
the spins parallel and so is subject to quantum fluctuations induced by the 
transverse part of the interaction. This physics will reappear when we study 
non-collinear states in QHE magnets away from fllling factor v = 1. 

The magnon dispersion for the ferromagnet can be understood in terms of 
bosonic 'particle' (the flipped spin) hopping on the lattice with a tight-binding 
model dispersion relation. The magnons are bosons because spin operators on 
different sites commute. They are not free bosons however because of the hard 
core constraint that (for spin 1/2) there can be no more than one flipped spin 
per site. Hence multi-magnon excited states can not be computed exactly. Some 
nice renormalization group arguments about magnon interactions can be found 
in [47]. 

The QHE ferromagnet is itinerant and we have to develop a somewhat differ- 
ent picture. Nevertheless there will be strong similarities to the lattice Heisen- 
berg model. The exact ground state is given by eq. (1.183) with 

$(^1, ...,zn)= Y[{z^ - zj) I-*"!'. (1.195) 

i<j 



^^We use the phase factor e "I'^j here rather than e"'"'^'^-' simply to be consistent with 



being the Fourier transform of S ■ 
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To find the spin wave excited states we need to find the analog of eq. (1.190). 
The Fourier transform of the spin lowering operator for the continuum system 
is 

N 

where fj is the position operator for the jth particle. Recall from eq. (1.144) 
that we had to modify Feynman's theory of the collective mode in superfluid 
helium by projecting the density operator onto the Hilbert space of the lowest 
Landau level. This suggests that we do the same in eq. (1.196) to obtain the 
projected spin flip operator. In contrast to the good but approximate result 
we obtained for the collective density mode, this procedure actually yields the 
exact one-magnon excited state (much like we found for the lattice model). 
Using the results of appendix A, the projected spin lowering operator is 

5- = e-^l«l^f^r,(j)5- (1.197) 

where q is the complex number representing the dimensionless wave vector qi 

and Tq{j) is the magnetic translation operator for the jth particle. The com- 
mutator of this operator with the Coulomb interaction Hamiltonian is 

[H,S-] = ^^t^(fc) [p-kPk,S-] 

fc#0 

= ^E^e^) {P-I^ [Pi^'^^] + ^4 • (1-198) 

We will shortly be applying this to the fully polarized ground state l'^). As 
discussed in appendix A, no density wave excitations are allowed in this state 
and so it is annihilated by pk- Hence we can without approximation drop the 
second term above and replace the first one by 

fc#0 

Evaluation of the double commutator following the rules in appendix A yields 
|V)=e,5- IV) (1.200) 

where 

eg = 2^e-il'=l%(A;) sm^(l-qAk). (1.201) 

Since is an eigenstate of H, this proves that S~ \tp) is an exact excited 
state of H with excitation energy Cg. In the presence of the Zeeman coupling 

Cq ^ eq + A. 



(1.196) 
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Figure 1.23: Schematic illustration of the QHE ferromagnet spinwave dispersion. 
There is a gap at small k equal to the Zeeman splitting, Az. At large wave 
vectors, the energy saturates at the Coulomb exchange energy scale A^; + Az ~ 
lOOK. 



This result tells us that, unlike the case of the density excitation, the single- 
mode approximation is exact for the case of the spin density excitation. The 
only assumption we made is that the ground state is fully polarized and has 
u=l. 

For small q the dispersion starts out quadratically 

eg ~ Aq^ (1.202) 

with 

A= ^^e-5l'=l%()fe) (1.203) 

as can be seen by expanding the sine function to lowest order. For very large q 
sin^ can be replaced by its average value of ^ to yield 

e,~^t;(fc)e- 5 (1.204) 

Thus the energy saturates at a constant value for g — !■ cx) as shown in fig. (1.23). 
(Note that in the lattice model the wave vectors are restricted to the first Bril- 
louin zone, but here they are not.) 

While the derivation of this exact result for the spin wave dispersion is alge- 
braically rather simple and looks quite similar (except for the LLL projection) 
to the result for the lattice Heisenberg model, it does not give a very clear 
physical picture of the nature of the spin wave collective mode. This we can 
obtain from eq. (1.197) by noting that Tq{j) translates the particle a distance 
q X zi"^ . Hence the spin wave operator S~ flips the spin of one of the partic;les 
and translates it spatially leaving a hole behind and creating a particle-hole 
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Figure 1.24: Illustration of the fact that the spin flip operator causes translations 
when projected into the lowest Landau level. For very large wave vectors the 
particles is translated completely away from the exchange hole and loses all its 
favorable Coulomb exchange energy. 

pair carrying net momentum proportional to their separation as illustrated in 
fig. (1.24). For large separations the cxcitonic Coulomb attraction between the 
particle and hole is negligible and the energy cost saturates at a value related to 
the Coulomb exchange energy of the ground state given in eq. (1.113). The ex- 
act dispersion relation can also be obtained by noting that scattering processes 
of the type illustrated by the dashed lines in fig. (1.24) mix together Landau 
gauge states 

I TTTTTT) (1.205) 

with different wave vectors k. Requiring that the state be an eigenvector of 
translation uniquely restricts the mixing to linear combinations of the form 

4-g„i %^ I TTTTTT). (1.206) 



Evaluation of the Coulomb matrix elements shows that this is indeed an exact 
eigenstate. 



1.10.4 Effective Action 

It is useful to try to reproduce these microscopic results for the spin wave ex- 
citations within an effective field theory for the spin degrees of freedom. Let 
m(r ) be a vector field obeying m-m = l which describes the local orientation of 
the order parameter (the magnetization). Because the Coulomb forces are spin 
independent, the potential energy cost can not depend on the orientation of to 
but only on its gradients. Hence we must have to leading order in a gradient 
expansion 

U =\ps j (fr di^m" d^m" - ^nA j cfr m' (1.207) 

where Ps is a phenomenological 'spin stiffness' which in two dimensions has units 
of energy and n = is the particle density. We will learn how to evaluate it 
later. 
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Wc can think of this expression for the energy as the leading terms in a 
functional Taylor series expansion. Symmetry requires that (except for the 
Zeeman term) the expression for the energy be invariant under uniform global 
rotations of m. In addition, in the absence of disorder, it must be translationally 
invariant. Clearly the expression in (1.207) satisfies these symmetries. The only 
zero-derivative term of the appropriate symmetry is m^m^ which is constrained 
to be unity everywhere. There exist terms with more derivatives but these 
are irrelevant to the physics at very long wavelengths. (Such terms have been 
discussed by Read and Sachdev [47].) 

To understand how time derivatives enter the effective action wc have to 
recall that spins obey a first-order (in time) precession equation under the in- 
fluence of the local exchange field. Consider as a toy model a single spin in 
an external field A. 

H = -nA"S" (1.208) 

The Lagrangian describing this toy model needs to contain a first order time 
derivative and so must have the form (see discussion in appendix B) 

jC = hS {-rh>'A^'[m] + A^m" + X{m^'m^' - 1)} (1.209) 

where 5 = | is the spin length and A is a Lagrange multiplier to enforce the 
fixed length constraint. The unknown vector A can be determined by requiring 
that it reproduce the correct precession equation of motion. The precession 
equation is 

= e"M/3^«5./3 (1.210) 

S = -AxS (1-211) 

which corresponds to counterclockwise precession around the magnetic field. 

We must obtain the same equation of motion from the Euler-Lagrange equa- 
tion for the Lagrangian in eq. (1.209) 

— - =0 (1.212) 

which may be written as 

A'' 2Am'' = F^'^'m'' (1.213) 

where 

F"- = d^A, - d^A^ (1.214) 



^^That is, the Coulomb exchange energy which tries to keep the spins locally parallel. In a 
Haxtree-Fock picture we could represent this by a term of the form —h{r) ■ s{r) where ) 
is the self-consistent field. 
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and 9p means {not the derivative with respect to some spatial coordinate). 
Since F'^" is antisymmetric let us guess a solution of the form 

F^"" = e">"'m'^. (1.215) 
Using this in eq. (1.213) yields 

A'' + 2Xm^' = e^^'m^m". (1.216) 
Applying e^^^m^ to both sides and using the identity 

,.a/3,.Ar, ^ ^^^^^^ _ ^^^^^^ (i_217) 

we obtain 

-{Kxrhy =m^ -rn'{m^m'^). (1.218) 

The last term on the right vanishes due to the length constraint. Thus we find 
that our ansatz in eq. (1.215) does indeed make the Euler-Lagrange equation 
correctly reproduce eq. (1.211). 
Eq. (1.215) is equivalent to 

Vm X = m (1.219) 

indicating that A is the vector potential of a unit magnetic monopole sitting 
at the center of the unit sphere on which m lives as illustrated in fig. (1.25). 
Note (the always confusing point) that we are interpreting rh as the coordinate 
of a fictitious particle living on the unit sphere (in spin space) surrounding the 
monopole. 

Recalling eq. (1.20), we see that the Lagrangian for a single spin in eq. (1.209) 

is equivalent to the Lagrangian of a massless object of charge —5, located at 
position rh, moving on the unit sphere containing a magnetic monopole. The 
Zeeman term represents a constant electric field —A producing a force AS on 
the particle. The Lorentz force caused by the monopole causes the particle to 
orbit the sphere at constant 'latitude'. Because no kinetic term of the form 
m"m" enters the Lagrangian, the charged particle is massless and so lies only 
in the lowest Landau level of the monopole field. Note the similarity here to 
the previous discussion of the high field limit and the semiclassical percolation 
picture of the integer Hall effect. For further details the reader is directed to 
appendix B and to Haldane's discussion of monopole spherical harmonics [48] . 

If the 'charge' moves slowly around a closed counterclockwise path m{t) 
during the time interval [0, T] as illustrated in fig. (1.25), the quantum amplitude 



j- dtC 

Jo 



(1.220) 



contains a Berry's phase [49] contribution proportional to the 'magnetic flux' 
enclosed by the path 

—iS f dtra^ A!^ —iS S A.-drri / 1 001 \ 

e Jo = e J . (1.221) 
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Figure 1.25: Magnetic monopole in spin space. Arrows indicate the curl of the 
Berry connection V x ^ emanating from the origin. Shaded region indicates 
closed path rn{t) taken by the spin order parameter during which it acquires 
a Berry phase proportional to the monopole flux passing through the shaded 
region. 



As discussed in appendix B, this is a purely geometric phase in the sense that 
it depends only on the geometry of the path and not the rate at which the path 
is traversed (since the expression is time reparameterization invariant). Using 
Stokes theorem and eq. (1.219) we can write the contour integral as a surface 
integral 

^-iS ^ A-dm _ ^-iS J dQ-VxA _ ^-iSQ ^-^ 222) 

where dQ = rhdfl is the directed area (solid angle) element and O is the to- 
tal solid angle subtended by the contour as viewed from the position of the 
monopole. Note from fig. (1.25) that there is an ambiguity on the sphere as to 
which is the inside and which is the outside of the contour. Since the total solid 
angle is 47r we could equally well have obtained^^ 

g+iS(47r-a)^ (1.223) 

Thus the phase is ambiguous unless 5* is an integer or half-integer. This consti- 
tutes a 'proof that the quantum spin length must be quantized. 

Having obtained the correct Lagrangian for our toy model we can now readily 
generalize it to the spin wave problem using the potential energy in eq. (1.207) 



C = -hSn J d^r I m^(f ) A''[m] - Am%r ] 



^^The change in the sign from -i-i to —i is due to the fact that the contour switches from 
being counterclockwise to clockwise if viewed as enclosing the 4n — CI area instead of the f2 
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-^Ps J d^r d^wyd^m" + j (fr A(f ) (m'^m" - 1). (1.224) 

The classical equation of motion can be analyzed just as for the toy model, 
however we will take a slightly different approach here. Let us look in the low 
energy sector where the spins all li(^ closc^ to the z direction. Then we can write 



m 



Now choose the 'symmetric gauge' 



which obeys eq. (1.219) for m close to z. 

Keeping only quadratic terms in the Lagrangian we obtain 

C = -nSn [ £r {Umym'' -rn'my) 



-A I 1 - -rrv^m 



1 



This can be conveniently rewritten by defining a complex field 



(1.225) 



(1.226) 



(1.227) 



-Snh 



r 



dt 



dt 



The classical equation of motion is the Schrodinger like equation 



dt nS ^ 
This has plane wave solutions with quantum energy 

nS 



(1.228) 



(1.229) 



(1.230) 



We can fit the phenomenological stiffness to the exact dispersion relation in 
eq. (1.202) to obtain 



nS 



(1.231) 



fe#0 




Figure 1.26: Illustration of a skyrmion spin texture. The spin is down at the 
origin and gradually turns up at infinite radius. At intermediate distances, the 
XY components of the spin exhibit a vortex-like winding. Unlike a C/(l) vortex, 
there is no singularity at the origin. 



Exercise 1.20 Derive eq. (1.231) from first principles by evaluating the loss 
of exchange energy when the Landau gauge v ^ I ground state is distorted to 
make the spin tumble in the x direction 




(1.232) 



where 6k = —^kl"^ and j = ^ is the (constant) spin rotation angle gradient 
(since x = —kP in this gauge). 



1.10.5 Topological Excitations 

So far we have studied neutral collective excitations that take the form of spin 
waves. They are neutral because as we have seen from eq. (1.197) they consist 
of a particle- hole pair. For very large momenta the spin- flipped particle is 
translated a large distance q x z(? away from its original position as discussed in 
appendix A. This looks locally like a charged excitation but it is very expensive 
because it loses all of its exchange energy. It is sensible to inquire if it is possible 
to make a cheaper charged excitation. This can indeed be done by taking into 
account the desire of the spins to be locally parallel and producing a smooth 
topological defect in the spin orientation [46, 50-56] known as a skyrmion by 
analogy with related objects in the Skyrme model of nuclear physics [57]. Such 
an object has the beautiful form exhibited in fig. (1.26). Rather than having 
a single spin suddenly flip over, this object gradually turns over the spins as 
the center is approached. At intermediate distances the spins have a vortex-like 
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configuration. However unlike a (7(1) vortex, there is no singularity in the core 
region because the spins are able to rotate downwards out of the xy plane. 

In nuclear physics the Skyrme model envisions that the vacuum is a 'fer- 
romagnct' described by a four component field $^ subject to the constraint 
^fi^fi — i_ There are three massless (i.e. linearly dispersing) spin wave exci- 
tations corresponding to the three directions of oscillation about the ordered 
direction. These three massless modes represent the three (nearly) massless pi- 
ous TT"*", 7r°, 7r~. The nucleons (proton and neutron) are represented by skyrmion 
spin textures. Remarkably, it can be shown (for an appropriate form of the ac- 
tion) that these objects are fermions despite the fact that they are in a sense 
made up of a coherent superposition of (an infinite number of) bosonic spin 
waves. 

We shall see a very similar phenomenology in QHE fcrromagnets. At filling 
factor u, skyrmions have charge ±i/e and fractional statistics much like Laughlin 
quasiparticles. For v = 1 these objects are fermions. Unlike Laughlin quasi- 
particles, skyrmions are extended objects, and they involve many flipped (and 
partially flipped) spins. This property has profound implications as we shall 
see. 

Let us begin our analysis by understanding how it is that spin tc^xtures can 
carry charge. It is clear from the Pauli principle that it is necessary to flip at 
least some spins to locally increase the charge density in a = 1 ferromagnet. 
What is the sufficient condition on the spin distortions in order to have a den- 
sity fluctuation? Remarkably it turns out to be possible, as we shall see, to 
uniquely express the charge density solely in terms of gradients of the local spin 
orientation. 

Consider a ferromagnet with local spin orientation m(r ) which is static. As 
each electron travels we assume that the strong exchange field keeps the spin 
following the local orientation m. If the electron has velocity x^, the rate of 
change of the local spin orientation it sees is m'^ = x^-^^m". This in turn 
induces an additional Berry's phase as the spin orientation varies. Thus the 
single-particle Lagrangian contains an additional first order time derivative in 
addition to the one induced by the magnetic field coupling to the orbital motion 

Co = --i'^A^ + hSm^Anm]. (1.233) 
c 

Here refers to the electromagnetic vector potential and A'^ refers to the 
monopole vector potential obeying eq. (1.219) and we have set the mass to zero 
(i.e. dropped the x^^x^^ term). This can be rewritten 

Cq = -^x^{A^' + a^') (1.234) 

where (with <I>o being the flux quantum) 

= -<^oS (^"i") -4"N (1-235) 

represents the 'Berry connection', an additional vector potential whic;li repro- 
duces the Berry phase. The additional fake magnetic flux due to the curl of the 



The Quantum Hall Effect 



73 



Berry connection is 



The first term vanishes by symmetry leaving 

where F"^ is given by eq. (1.215) and we have taken advantage of the fact that 
the remaining factors are antisymmetric under the exchange z/ -(-^ 7. Using 
eq. (1.215) and setting 5 = ^ we obtain 

b = -^oP (1.238) 

where 

p = - — e '^e m CqTO Ogm 

= -^^"'^m ■ 5„m X d/sm (1.239) 

is (for reasons that will become clear shortly) called the topological density or 
the Pontryagin density. 

Imagine now that we adiabatically deform the uniformly magnetized spin 
state into some spin texture state. We see from eq. (1.238) that the orbital 
degrees of freedom see this as adiabatically adding additional flux b{f). Recall 
from eq. (1.171) and the discussion of the charge of the Laughlin quasiparticle, 
that extra charge density is associated with extra flux in the amount 

dp = ^a,yb (1.240) 
Sp = vep. (1.241) 

Thus we have the remarkable result that the changes in the electron charge 
density are proportional to the topological density. 

Our assumption of adiabaticdty iw valid as long as the spin fluctuation fre- 
quency is much lower than the charge excitation gap. This is an excellent 
approximation for = 1 and still good on the stronger fractional Hall plateaus. 

It is interesting that the fermionic charge density in this model can be ex- 
pressed solely in terms of the vector boson field rn{f), but there is something 
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even more significant here. The skyrmion spin texture has total topological 

charge 

Qtop =^J(fr e"^m ■ d^m x d^m (1.242) 

which is always an integer. In fact for any smooth spin texture in which the 
spins at infinity are all parallel, Qtop is always an integer. Since it is impossible 
to continuously deform one integer into another, Qtop is a topological invariant. 
That is, if Qtop = ±1 because a skyrmion (anti-skyrmion) is present, Qtop is 
stable against smooth continuous distortions of the field rri. For example a 
spin wave could pass through the skyrmion and Qtop would remain invariant. 
Thus this charged object is topologically stable and has fermion number (i.e., 
the number of fermions (electrons) that flow into the region when the object is 
formed) 

N = pQtop. (1.243) 

For = 1, A'^ is an integer (±1 say) and has the fermion number of an electron. 
It is thus continuously connected to the single flipped spin example discussed 
earlier. 

We are thus led to the remarkable conclusion that the spin degree of freedom 
couples to the electrostatic potential. Because skyrmions carry charge, we can 
aff'ect the spin configuration using electric rather than magnetic fields! 

To understand how Qtop always turns out to be an integer, it is useful to 
consider a simpler case of a one-dimensional ring. We follow here the discussion 
of [58]. Consider the unit circle (known to topologists as the one-dimensional 
sphere Si). Let the angle 9 e[0, 27r] parameterize the position along the curve. 
Consider a continuous, suitably well-behaved, complex function ip^O) = e*'^^^^ 
defined at each point on the circle and obeying = 1. Thus associated with 
each point 9 is another tmit circle giving the possible range of values of ip{9). 
The function ip{9) thus defines a trajectory on the torus x 5*1 illustrated in 
fig. (1.27). The possible functions ip{6) can be classifled into diff'erent homotopy 
classes according to their winding number n G Z 

Because the points 9 = and 9 = 2tt are identified as the same point 

V'(O) = -0(277) ^ ^(27r) - 1^(0) = 27r X integer (1.245) 

and so n is an integer. Notice the crucial role played by the fact that the 
'topological density' ^ ^ is the Jacobian for converting from the coordinate 
d in the domain to the coordinate ip in the range. It is this fact that makes the 
integral in eq. (1.244) independent of the detailed local form of the mapping 
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Figure 1.27: Illustration of mappings (fi{6) with: zero winding number (left) and 
winding number +2 (right). 



(p{0) and depend only on the overall winding number. As we shall shortly see, 
this same feature will also turn out to be true for the Pontryagin density. 

Think of the function (p{9) as defining the path of an elastic band wrapped 
around the torus. Clearly the band can be stretched, pulled and distorted in 
any smooth way without any effect on n. The only way to change the winding 
number from one integer to another is to discontinuously break the elastic band, 
unwind (or wind) some extra turns, and then rejoin the cut pieces. 

Another way to visualize the homotopy properties of mappings from Si to 
5*1 is illustrated in fig. (1.28). The solid circle represents the domain and the 
dashed circle represents the range ip. It is useful to imagine the 6 circle as being 
an elastic band (with points on it labeled by coordinates running from to 2n) 
which can be 'lifted up' to the (p circle in such a way that each point of 6 lies just 
outside the image point ip{9). The figure illustrates how the winding number n 
can be interpreted as the number of times the domain 9 circle wraps around the 
range ip circle. (Note: even though the elastic band is 'stretched' and may wrap 
around the ip circle more than once, its coordinate labels still only run from 
to 27r.) This interpretation is the one which we will generalize for the case of 
skyrmions in 2D ferromagnets. 

We can think of the equivalence class of mappings having a given winding 
number as an element of a group called the homotopy group 7ri(5'i). The group 
operation is addition and the winding number of the sum of two functions, 
ip{9) = ipi{9) + ip2{9), is the sum of the two winding numbers n = ni+n2. Thus 
7ri(S'i) is isomorphic to Z, the group of integers under addition. 

Returning now to the ferromagnet we see that the unit vector order parame- 
ter m defines a mapping from the plane i?2 to the two-sphere 5*2 (i.e. an ordinary 
sphere in three dimensions having a two-dimensional surface). Because we as- 
sume that m = z for all spatial points far from the location of the skyrmion, 
we can safely use a projective map to 'compactify' R2 into a sphere 5*2. In 
this process all points at infinity in R2 are mapped into a single point on ^2, 
but since fh{f) is the same for all these different points, no harm is done. We 
are thus interested in the generalization of the concept of the winding number 
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Figure 1.28: A different representation of the mappings from 6 to (p. The 
dashed line represents the domain 9 and the sohd hne represents the range (p. 
The domain is 'hfted up' by the mapping and placed on the range. The winding 
number n is the number of times the dashed circle wraps the solid circle (with 
a possible minus sign depending on the orientation). 
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Figure 1.29: Infinitesimal circuit in spin space associated with an infinitesimal 
circuit in real space via the mapping m(r). 



to the mapping ^2 82- The corresponding homotopy group 712(32) is also 
equivalent to Z as we shall see. 

Consider the following four points in the plane and their images (illustrated 
in fig. (1.29)) under the mapping 



ix,y) 

{x + dx, y) 
[x, y + dy) 
{x + dx,y + dy) 



m{x, y) 
m{x + dx, y) 
m{x, y + dy) 
rh{x + dx, y + dy). 



(1.246) 



The four points in the plane define a rectangle of area dxdy. The four points on 
the order parameter (spin) sphere define an approximate parallelogram whose 
area (solid angle) is 



dw 



[m{x + dx, y) - rh{x, y)] x [m{x, y + dy) - rh{x, y)] ■ m{x, y) 



m ■ dufh X dvfh dxdy 



A-Kp dxdy. 



(1.247) 



Thus the Jacobian converting area in the plane into solid angle on the sphere 
is 47r times the Pontryagin density p. This means that the total topological 
charge given in eq. (1.242) must be an integer since it counts the number of 
times the compactified plane is wrapped around the order parameter sphere by 
the mapping. The 'wrapping' is done by lifting each point r in the compactified 
plane up to the corresponding point rh{r) on the sphere just as was described 
for 7ri(S'i) in fig. (1.28). 

For the skyrmion illustrated in fig. (1.26) the order parameter function m(r ) 
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was chosen to be the standard form that minimizes the gradient energy [58] 

= (1.248a) 

my = ^^^^^i';^) (1.248b) 
r^-A^ 

= (1.248c) 

where (r, 0) are the polar coordinates in the plane, A is a constant that controls 
the size scale, and <^ is a constant that controls the XY spin orientation. (Ro- 
tations about the Zeeman axis leave the energy invariant.) Prom the figure it is 
not hard to see that the skyrmion mapping wraps the compactified plane around 
the order parameter sphere exactly once. The sense is such that Qtop = — 1- 



Exercise 1.21 Show that the topological density can be written in polar spa- 
tial coordinates as 

1 _ din dm 
^ 47rr"^ dr 86 
Use this result to show 

1 / 2A 



47rVA2+r2 
and hence 

Qtop = — 1 

for Ui( skijnnion viapphuj rii cqs. (l.i'JiHa l.'iJi.Hc). 



It is worthwhile to note that it is possible to write down simple microscopic 

variational wave functions for the skyrmion which arc closely related to the 
continuum field theory results obtained above. Consider the following state in 
the plane [51] 



A 



*i, (1.249) 



where is the v = 1 filled Landau level state refers to the spinor for 
the j'th particle, and A is a fixed length scale. This is a skyrmion because it 
has its spin purely down at the origin (where Zj = 0) and has spin purely up 
at infinity (where \zj\ 3> A). The parameter A is simply the size scale of the 
skyrmion [46,58]. At radius A the spinor has equal weight for up and down spin 
states (since \zj\ = A) and hence the spin lies in the XY plane just as it does 
for the solution in eq. (1.248c). Notice that in the limit A — > (where the 
continuum eflfective action is invalid but this microscopic wave function is still 
sensible) we recover a fully spin polarized filled Landau level with a charge- 1 
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Laughlin quasiholc at the origin. Hence the number of flipped spins interpolates 
continuously from zero to infinity as A increases. 

In order to analyze the skyrmion wave function in eq. (1.249), we use the 
Laughlin plasma analogy. Recall from our discussion in sec. 1.6.1 that in this 
analogy the norm of tp\, Tr^^j J D[z] \^[z]\^ is viewed as the partition function 
of a Coulomb gas. In order to compute the density distribution we simply need 
to take a trace over the spin 



This partition function describes the usual logarithmically interacting Coulomb 
gas with uniform background charge plus a spatially varying impurity back- 
ground charge Apb{r), 



AP.M - -^VV(0=+^^(^4^, (1.251) 
V{r) = -llog(r2 + A2). (1.252) 



For large enough scale size X ^ £, local neutrality of the plasma [59] forces 
the electrons to be expelled from the vicinity of the origin and implies that 
the excess electron number density is precisely — Apb(r), so that eq. (1.251) is 
in agreement with the standard result [58] for the topological density given in 
ex. 1.21. 

Just as it was easy to find an explicit wave function for the Laughlin quasi- 
hole but proved difficult to write down an analytic wave function for the Laughlin 
quasi-electron, it is similarly difficult to make an explicit wave function for the 
anti-skyrmion. Finally, we note that by replacing (^) by (^„) , we can generate 
a skyrmion with a Pontryagin index n. 
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Exercise 1.22 The argument given above for the charge density of the mi- 
croscopic skyrmion state wave function used local neutrality of the plasma and 
hence is valid only on large length scales and thus requires X ^ i. Find the 
complete microscopic analytic solution for the charge density valid for arbi- 
trary X, by using the fact that the proposed manybody wave function is nothing 
but a Slater determinant of the single particle states (j)m{z), 

4>m{z) = ('^ e-^ . (1.253) 



^27r2™+im! (m + 1 + ^) 



Show that the excess electron number density is then 

AnW(^)^ (1-254) 

m=0 

which yields 

An^^\z) = i- Q da e--^^^-°'\\z\'' + X^) - 1^ . (1.255) 

Similarly, find the spin density distribution (r) and show that it also agrees 
with the field-theoretic expression in eq. (1.248c) in the large X limit. 



The skyrmion solution in eqs. (1.248a-1.248c) minimizes the gradient energy 
Eo = \ps J d\ d^nfd^m". (1.256) 

Notice that the energy cost is scale invariant since this expression contains two 

integrals and two derivatives. Hence the total gradient energy is independent of 
the scale factor A and for a single skyrmion is given by [46, 58] 

^0 = 47rp« = ^e<x, (1.257) 

where e<x> is the asymptotic large q limit of the spin wave energy in eq. (1.201). 

Since this spin wave excitation produces a widely separated particle-hole pair, we 
see that the energy of a widely separated skyrmion- antiskyrmion pair ( j + 2) ^oo 
is only half as large. Thus skyrmions are considerably cheaper to create than 
simple flipped spins. 

Notice that eq. (1.257) tells us that the charge excitation gap, while only 
half as large as naively expected, is finite as long as the spin stiffness ps is finite. 
Thus we can expect a dissipationless Hall plateau. Therefore, as emphasized by 
Sondhi et al. [46] , the Coulomb interaction plays a central role in the v = 1 inte- 
ger Hall effect. Without the Coulomb interaction the charge gap would simply 
be the tiny Zeeman gap. With the Coulomb interaction the gap is large even 

■^''This energy advantage is reduced if the finite thickness of the inversion layer is taken into 
account. The skyrmion may in some cases turn out to be disadvantageous in higher Landau 
levels. 
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in the limit of zero Zccman energy because of the spontaneous ferromagnetic 
order induced by the spin stiffness. 

At precisely u = 1 skyrmion/antiskyrmion pairs will be thermally activated 
and hence exponentially rare at low temperatures. On the other hand, because 
they are the cheapest way to inject charge into the system, there will be a 
finite density of skyrmions even in the ground state if v ^ 1. Skyrmions also 
occur in ordinary 2D magnetic films but since they do not carry charge (and 
are energetically expensive since ps is quite large) they readily freeze out and 
are not particularly important. 

The charge of a skyrmion is sharply quantized but its number of flipped 
spins depends on its area ~ A^. Hence if the energy were truly scale invariant, 
the number of flipped spins could take on any value. Indeed one of the early 
theoretical motivations for skyrmions was the discovery in numerical work by 
Rezayi [46,60] that adding a single charge to a filled Landau level converted the 
maximally ferromagnetic state into a spin singlet. In the presence of a finite 
Zeeman energy the scale invariance is lost and there is a term in the energy that 
scales with AA^ and tries to minimize the size of the skyrmion. Competing with 
this however is a Coulomb term which we now discuss. 

The Lagrangian in eq. (1.224) contains the correct leading order terms in 
a gradient expansion. There are several possible terms which are fourth order 
in gradients, but a particular one dominates over the others at long distances. 
This is the Hartree energy associated with the charge density of the skyrmion 



and e is the dielectric constant. The long range of the Coulomb interaction 
makes this effectively a three gradient term that distinguishes it from the other 
possible terms at this order. Recall that the Coulomb interaction already entered 
in lower order in the computation of ps- That however was the exchange energy 
while the present term is the Hartree energy. The Hartree energy scales like 

2 

and so prefers to expand the skyrmion size. The competition between the 
Coulomb and Zeeman energies yields an optimal number of approximately four 
flipped spins according to microscopic Hartree Fock calculations [61]. 

Thus a significant prediction for this model is that each charge added (or 
removed) from a filled Landau level will flip several (~ 4) spins. This is very 
different from what is expected for non-interacting electrons. As illustrated in 
fig. (1.30) removing an electron leaves the non-interacting system still polarized. 
The Pauli principle forces an added electron to be spin reversed and the mag- 
netization drops from unity at = 1 to zero at = 2 where both spin states of 
the lowest Landau level are fully occupied. 

Direct experimental evidence for the existence of skyrmions was first ob- 
tained by Barrett et al. [62] using a novel optically pumped NMR technique. 




(1.258) 



where 



5p 




(1.259) 
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I 

Figure 1.30: Illustration of the spin configurations for non-interacting electrons 
at filling factor v = 1 in the presence of a hole (top) and an extra electron 
(bottom). 



The Hamiltonian for a nucleus is [63] 

Hn = -AnP + 07 • s (1.260) 

where / is the nuclear angular momentum. A^v is the nuclear Zccman frequency 
(about 3 orders of magnitude smaller than the electron Zeeman frequency), O 
is the hyperfine coupling and s is the electron spin density at the nuclear site. 
If, as a first approximation we replace s by its average value 

Hn ^ {-An + ^(s^)) r (1.261) 

we see that the precession frequency of the nucleus will be shifted by an amount 
proportional to the magnetization of the electron gas. The magnetization de- 
duced using this so-called Knight shift is shown in fig. (1.31). The electron 
gas is 100% polarized a,t v = 1, but the polarization drops off sharply (and 
symmetrically) as charge is added or subtracted. This is in sharp disagreement 
with the prediction of the free electron model as shown in the figure. The ini- 
tial steep slope of the data allows one to deduce that 3.5-4 spins reverse for 
each charge added or removed. This is in excellent quantitative agreement with 
Hartrce-Fock calculations for the skyrmion model [61]. 

Other evidence for skyrmions comes from the large change in Zeeman energy 
with field due to the large number of flipped spins. This has been observed in 
transport [64] and in optical spectroscopy [65]. Recall that spin-orbit effects in 
GaAs make the electron g factor —0.4. Under hydrostatic pressure g can be 
tuned towards zero which should greatly enhance the skyrmion size. Evidence 
for this effect has been seen [66]. 
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Figure 1.31: NMR Knight shift measurement of the electron spin polarization 
near filling factor u = I. Circles arc the data of Barrett et al. [62]. The dashed 
line is a guide to the eye. The solid line is the prediction for non-interacting 
electrons. The peak represents 100% polarization at v = 1. The steep slope on 
each side indicates that many 4) spins flip over for each charge added (or 
subtracted). The observed symmetry around f = 1 is due to the particle- hole 
symmetry between skyrmions and antiskyrmions not present in the free-electron 
model. 
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Figure 1.32: NMR nuclear spin relaxation rate 1/Ti as a function of filling 

factor. After Tycko et al. [68]. The relaxation rate is very small at = 1, but 
rises dramatically away from v = 1 due to the presence of skyrmions. 

1.11 Skyrmion Dynamics 

NMR [62] and nuclear specific heat [67] data indicate that skyrmions dramat- 
ically enhance the rate at which the nucl(;ar spins relax. This nuclear spin 
relaxation is due to the transverse terms in the hyperfine interaction which we 
neglected in discussing the Knight shift 

Ul{I+s-+I-s+) = hl |7+5^5r+h.c.|. (1.262) 

The free electron model would predict that it would be impossible for an electron 
and a nucleus to undergo mutual spin flips because the Zeeman energy would not 
be conserved. (Recall that A^v ~ The spin wave model shows that the 

problem is even worse than this. Recall from eq. (1.201) that the spin Coulomb 
interaction makes spin wave energy much larger than the electron Zeeman gap 
except at very long wavelengths. The lowest frequency spin wave excitations 
lie above 20-50 GHz while the nuclei precess at 10-100 MHz. Hence the two 
sets of spins are unable to couple effectively. Ki v = 1 this simple picture is 
correct. The nuclear relaxation time Ti is extremely long (tens of minutes to 
many hours depending on the temperature) as shown in fig. (1.32). However 
the figure also shows that iov v ^ 1 the relaxation rate 1/Ti rises dramatically 
and Ti falls to ~ 20 seconds. In order to understand this dramatic variation we 
need to develop a theory of spin dynamics in the presence of skyrmions. 
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The 1/Ti data is telling us that for 7^ 1 at least some of the electron 
spin fluctuations are orders of magnitude lower in frequency than the Zeeman 
splitting and these low frequency modes can couple strongly to the nuclei. One 
way this might occur is through the presence of disorder. We see from cq. (1.262) 
that NMR is a local probe which couples to spin flip excitations at all wave 
vectors. Recall from eq. (1.197) that lowest Landau level projection implies that 

contains a translation operator Tq. In the presence of strong disorder the 
Zeeman and exchange cost of the spin flips could be compensated by translation 
to a region of lower potential energy. Such a mechanism was studied in [69] but 
does not show sharp features in 1/Ti around v = I. 

We are left only with the possibility that the dynamics of skyrmions somehow 
involves low frequency spin fluctuations. For simplicity we will analyze this 
possibility ignoring the effects of disorder, although this may not be a valid 
approximation. 

Let us begin by considering a ferromagnetic u = 1 state containing a single 
skyrmion of the form parameterized in eqs. (1.248a 1.248c). There are two de- 
generacies at the classical level in the effective field theory: The energy does not 
depend on the position of the skyrmion and it does not depend on the angu- 
lar orientation ip. These continuous degeneracies are known as zero modes [58] 
and require special treatment of the quantum fluctuations about the classical 
solution. 

In the presence of one or more skyrmions, the quantum Hall ferromagnet 
is non-colinear. In an ordinary ferromagnet where all the spins are parallel, 
global rotations about the magnetization axis only change the quantum phase 
of the state — they do not produce a new state. Because the skyrmion has 
distinguishable orientation, each one induces a new U{1) degree of freedom in 
the system. In addition because the skyrmion has a distinguishable location, 
each one induces a new translation degree of freedom. As noted above, both 
of these are zero energy modes at the classical level suggesting that they might 
well be the source of low energy excitations which couple so effectively to the 
nuclei. We shall see that this is indeed the case, although the story is somewhat 
complicated by the necessity of correctly quantizing these modes. 

Let us begin by flnding the effective Lagrangian for the translation mode [8]. 
We take the spin conflguration to be 

m(f, t) = mo (f- R{t)^ (1.263) 

where mo is the static classical skyrmion solution and R{t) is the position degree 
of freedom. We ignore all other spin wave degrees of freedom since they arc 
gapped. (The gapless U{1) rotation mode will be treated separately below.) 
Eq. (1.224) yields a Berry phase term 

£0 = -hS j cPr m'^A^'lrn ] n{r ) (1.264) 

Rotation about the Zeeman alignment axis is accomplished hy R = e~'^'*'^ . But a 
colinear ferromagnet ground state is an eigenstate of S^, so rotation leaves the state invariant 
up to a phase. 
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where 

m'' = -^^^m^(f-^) (1.265) 

and unlike in eq. (1.224) we have taken into account our new-found knowledge 
that the density is non-uniform 

n{f) = no + ^e^"" m- df,rn X d^m. (1.266) 

OTT 

The second term in eq. (1.266) can be shown to produce an extra Berry phase 
when two skyrmions are exchanged leading to the correct minus sign for Fermi 
statistics (on the v = 1 plateau) but we will not treat it further. Eq. (1.264) 
then becomes 

Co = +nR''a''{f) (1.267) 
where the 'vector potential' 

o^(r ) = Sno j (fr (a^m")^" (1.268) 

has curl 



-Sno e^" J dPr dx {(O^m")^"} 

-Sno e^"" j dPr {d^m^) {dxm^) 



dm-y 



2 

= -27rnoQtop (1-269) 

Thus eq. (1.267) corresponds to the kinetic Lagrangian for a massless particle 
of charge — eQtop moving in a uniform magnetic field of strength B = But 
this of course is precisely what the skyrmion is [8]. 

We have kept here only the lowest order adiabatic time derivative term in 
the action. This is justific'd by the; existence of the spin excitation gap and 
the fact that we are interested only in much lower frequencies (for the NMR). 

If we ignore the disorder potential then the kinetic Lagrangian simply leads 
to a Hamiltonian that yields quantum states in the lowest Landau level, all 
of which are degenerate in energy and therefore capable of relaxing the nu- 
clei (whose precession frequency is extremely low on the scale of the electronic 
Zeeman energy). 

Let us turn now to the rotational degree of freedom represented by the co- 
ordinate (p in eqs. (1.248a-1.248c). The full Lagrangian is complicated and 

There may exist higher-order time-derivative terms which give the skyrmion a mass and 
there will also be damping due to radiation of spin waves at higher velocities. [70] 
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contains the degrees of freedom of the continuous field m{'r). We need to intro- 
duce the collective coordinate describing the orientation of the skyrmion as 
one of the degrees of freedom and then carry out the Feynman path integration 
over the quantum fluctuations in all the infinite number of remaining degrees of 
freedom. ^'^ This is a no n- trivial task, but fortunately we do not actually have to 
carry it out. Instead we will simply write down the answer. The answer is some 
functional of the path for the single variable (p{t). We will express this functional 
(using a functional Taylor series expansion) in the most general form possible 
that is consistent with the symmetries in the problem. Then we will attempt to 
identify the meaning of the various terms in the expansion and evaluate their 
coefficients (or assign them values phenomenologically) . After integrating out 
the high frequency spin wave fluctuations, the lowest-order symmetry-allowed 
terms in the action are 

C^ = hK<p+—p^ + ... (1.270) 

Again, there is a first-order term allowed by the lack of time-reversal symmetry 
and we have included the leading non-adiabatic correction. The full action 
involving m(r, t) contains only a first-order time derivative but a second order 
term is allowed by symmetry to be generated upon integrating out the high 
frequency fluctuations. We will not perform this explicitly but rather treat U 
as a phenomenological fitting parameter. 

The coefficient K can be computed exactly since it is simply the Berry phase 
term. Under a slow rotation of all the spins through 27r the Berry phase is (using 
eq. (B.22) in appendix B) 

J d^r n(r ) (-6'27r) [1 - mg(r )] = ^ ^ = '^'^K. (1.271) 

(The non-adiabatic term gives a 1/T contribution that vanishes in the adiabatic 
limit T ^ 00.) Thus we arrive at the important conclusion that K is the 
expectation value of the number of overturned spins for the classical solution 
mo{f). We emphasize that this is the Hartree-Fock (i.e., 'classical') skyrmion 
solution and therefore K need not be an integer. 

The canonical angular momentum conjugate to </? in eq. (1.270) is 

L,= ^-^=nK+^ip (1.272) 

and hence the Hamiltonian is 

Hip = L^ip — £<p 



■^^Examples of how to do this are discussed in various field theory texts, including Rajaraman 
[58]. 
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Having identified the Hamiltonian and expressed it in terms of the eoordinate 
and the canonical momentum conjugate to that coordinate, we quantize by 
simply making the substitution 

-ih-^ (1.274) 

to obtain 

= A (1.275) 

This can be interpreted as the Hamiltonian of a (charged) XY quantum rotor 
with moment of inertia /U circling a solenoid containing K flux quanta. (The 
Berry phase term in eq. (1.270) is then interpreted as the Aharonov-Bohm 
phase.) The eigenfunctions are 

V„(v) = ^e^™*' (1.276) 

V ZTT 



and the eigenvalues are 

em = ^{m-Kf. (1.277) 

The angular momentum operator is actually the operator giving the number 

of flipped spins in the skyrmion. Because of the rotational symmetry about the 
Zeeman axis, this is a good quantum number and therefore takes on integer val- 
ues (as required in any quantum system of finite size with rotational symmetry 
about the z axis) . The ground state value of m is the nearest integer to K. The 
ground state angular velocity is 

Hence if K is not an integer the skyrmion is spinning around at a finite veloc- 
ity. In any case the actual orientation angle ip for the skyrmion is completely 
uncertain since from eq. (1.276) 

l^Mf = ^ (1-279) 

(fi has a flat probability distribution (due to quantum zero point motion). We 
interpret this as telling us that the global U(l) rotation symmetry broken in 
the classical solution is restored in the quantum solution because of quantum 
fluctuations in the coordinate ip. This issue will arise again in our study of the 
Skyrme lattice where we will find that for an infinite array of skyrmions, the 
symmetry can sometimes remain broken. 

Microscopic analytical [71] and numerical [61] calculations do indeed find 
a family of low energy excitations with an approximately parabolic relation 
between the energy and the number of fiipped spins just as is predicted by 
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Figure 1.33: Electronic structure of the skyrniion lattice as determined by nu- 
merical Hartree-Fock calculations for filling factor v = 1.1 and Zeeman en- 
ergy 0.015|^. (a) Excess charge density (in units of l/(27r£^)) and (b) Two- 
dimensional vector representation of the XY components of the spin density. 
The spin stiffness makes the square lattice more stable than the triangular lat- 
tice at this filling factor and Zeeman coupling. Because of the U{1) rotational 
symmetry about the Zeeman axis, this is simply one representative member of 
a continuous family of degenerate Hartree-Fock solutions. After Brey et al. [71]. 

eq. (1.277). As mentioned earlier, ~ 4 for typical parameters. Except for 
the special case where K is a half integer the spectrum is non-degenerate and 
has an excitation gap on the scale of U which is in turn some fraction of the 
Coulomb energy scale ~ 100 K. In the absence of disorder even a gap of only 
1 K would make these excitations irrelevant to the NMR. We shall see however 
that this conclusion is dramatically altered in the case where many skyrmions 
are present. 

1.11.1 Skyrme Lattices 

For filling factors slightly away from v = 1 there will be a finite density of 
skyrmions or antiskyrmions (all with the same sign of topological charge) in the 
ground state [56,72.73]. Hartree-Fock calculations [72] indicate that the ground 
state is a Skyrme crystal. Because the skyrmions are charged, the Coulomb 
potential in eq. (1.258) is optimized for the triangular lattice. This is indeed the 
preferred structure for very small values of ]f — 1] where the skyrmion density is 
low. However at moderate densities the square lattice is preferred. The Hartree- 
Fock ground state has the angular variable (fij shifted by n between neighboring 
skyrmions as illustrated in fig. (1.33). This 'antiferromagnetic' arrangement 
of the XY spin orientation minimizes the spin gradient energy and would be 
frustrated on the triangular lattice. Hence it is the spin stiffness that stabilizes 
the square lattice structure. 

The Hartree-Fock ground state breaks both global translation and global 
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U{1) spin rotation symmetry. It is a kind of 'supersolid' with both diagonal 

= (s^(f ) s'{r')) (1.280) 
and off-diagonal 

G^ = {s+{f) s-ir')) (1.281) 

long-range order. For the case of a single skyrmion we found that the C/(l) 
symmetry was broken at the Hartree-Fock (classical) level but fully restored by 
quantum fluctuations of the zero mode coordinate ip. In the thermodynamic 
limit of an infinite number of skyrmions coupled together, it is possible for the 
global U{1) rotational symmetry breaking to survive quantum fluctuations.^^ 
If this occurs then an excitation gap is not produced. Instead we have a new 
kind of gapless spin wave Goldstone mode [74,75]. This mode is gapless despite 
the presence of the Zeeman field and hence has a profound effect on the NMR 
relaxation rate. The gapless Goldstone mode associated with the broken trans- 
lation symmetry is the ordinary magneto-phonon of the Wigner crystal. This 
too contributes to the nuclear relaxation rate. 

In actual practice, disorder will be important. In addition, the NMR ex- 
periments have so far been performed at temperatures which are likely well 
above the lattice melting temperature. Nevertheless the zero temperature lat- 
tice calculations to be discussed below probably capture the essential physics 
of this non co-linear magnet. Namely, there exist spin fluctuations at frequen- 
cies orders of magnitude below the Zeeman gap. At zero temperature these are 
coherent Goldstone modes. Above the lattice melting temperature they will be 
overdamped diffusive modes derived from the Goldstone modes. The essential 
physics will still be that the spin fluctuations have strong spectral density at 
frequencies far below the Zeeman gap. 

It turns out that at long wavelengths the magnetophonon and U{1) spin 
modes are decoupled. We will therefore ignore the positional degrees of freedom 
when analyzing the new U{1) mode. We have already found the U{1) Hamil- 
tonian for a single skyrmion in eq. (1.275). The simplest generalization to the 
Skyrme lattice which is consistent with the symmetries of the problem is 

H=^Y.^K,- Kf - cos {ipi - (1.282) 

where Kj = —ij^ is the angular momentum operator. The global U (1) symme- 
try requires that the interactive term be invariant if all of the ipj's are increased 
by a constant. In addition H must be invariant under cpj ipj + 2t: for any 
single skyrmion. We have assumed the simplest possible near-neighbor cou- 
pling, neglecting the possibility of longer range higher-order couplings of the 
form cosn(i^i — (pj) which are also symmetry allowcid. The phenomenological 

^^Loosely speaking this corresponds to the infinite system having an infinite moment of 
inertia (for global rotations) which allows a quantum wave packet which is initially localized 
at a pajrticulaj: orientation (f not to spread out even for long times. 
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coupling J must be negative to be consistent with the 'antifcrroniagnctic' XY 
order found in the Hartree-Fock ground state iUustrated in fig. (1.33). However 
we will find it convenient to instead make J positive and compensate for this 
by a 'gauge' change (pj ipj + tt on one sublattice. This is convenient because 
it makes the coupling 'ferromagnetic' rather than 'antiferromagnetic' 

Eq. (1.282) is the Hamiltonian for the quantum XY rotor model, closely 
related to the boson Hubbard model [76-78]. Readers familiar with super- 
conductivity will recognize that this model is commonly used to describe the 
superconductor- insulator transition in Josephson arrays [76,77]. The angular 
momentum eigenvalue of the Kj operator represents the number of bosons 
(Cooper pairs) on site j and the U term describes the charging energy cost 
when this number deviates from the electrostatically optimal value of K. The 
boson number is non-negative while Kj has negative eigenvalues. However we 
assume that K ^ 1 so that the negative angular momentum states are very 
high in energy. 

The J term in the quantum rotor model is a mutual torque that transfers 
units of angular momentum between neighboring sites. In the boson language 
the wave function for the state with m bosons on site j contains a factor 

V'm('^j) = e™*'^. (1.283) 

The raising and lowering operators are thus^^ 6^*"^^ . This shows us that the 
cosine term in eq. (1.282) represents the Josephson coupling that hops bosons 
between neighboring sites. 

For U ^ J the system is in an insulating phase well-described by the wave 
function 

V;(^i,(p2,...,<Pw) = ne""''' (1-284) 
j 

where m is the nearest integer to K. In this state every rotor has the same fixed 
angular momentum and thus every site has the same fixed particle number in 
the boson language. There is a large excitation gap 

Af^U{l-2\m- K\) (1.285) 

and the system is insulating. 

Clearly [■^l^ « 1 in this phase and it is therefore quantum disordered. That 
is, the phases {(fj} are wildly fluctuating because every configuration is equally 
likely. The phase fluctuations are nearly uncorrelated 

^giyj g-ivfc^ ^ g_|f^._f^^|/5_ (1.286) 



■^^ These operators have matrix elements (tpm+i\e~^'^'^\i)m) = 1 whereas a boson raising 
operator would have matrix element v^mT+T. For K ^ 1, m K and this is nearly a constant. 
Arguments like this strongly suggest that the boson Hubbard model and the quantum rotor 
model are essentially equivalent. In particular their order/disorder transitions are believed to 
be in the same universality class. 

^^An exception occurs if \m — K\ = ^ where the gap vanishes. See [78]. 
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For J ^ U the phases on neighboring sites are strongly coupled together 
and the system is a superconductor. A crude variational wave function that 
captures the essential physics is 

^{^„^„...,^^)^e'^<^r,-^^^-^^^ (1.287) 

where A is a variational parameter [79]. This is the simplest ansatz consistent 
with invariance under (pj ipj + 27r. For J ^ U, 1 and is large only 
for spin configurations with all of the XY spins locally parallel. Expanding the 
cosine term in eq. (1.282) to second order gives a harmonic Hamiltonian which 
can be exactly solved. The resulting gapless 'spin waves' are the Goldstone 
modes of the superconducting phase. 

For simplicity we work with the Lagrangian rather than the Hamiltonian 

r 1 

t^ = Y. ^^'^J + ^^^j +-/$^cos(<^i-^j) (1.288) 

3 (ij) 



The Berry phase term is a total derivative and can not affect the equations 
of motion.^^ Dropping this term and expanding the cosine in the harmonic 
approximation yields 

i (ij) 

This 'phonon' model has linearly dispersing gapless collective modes at small 
wavevectors 

huiq = VUJ qa (1.290) 

where a is the lattice constant. The parameters U and J can be fixed by 
fitting to microscopic Hartree-Fock calculations of the spin wave velocity and 
the magnetic susceptibility ('boson compressibility') [61,75]. This in turn allows 
one to estimate the regime of filling factor and Zeeman energy in which the U{1) 
symmetry is not destroyed by quantum fiuctuations [75] . 

Let us now translate all of this into the language of our non-colinear QHE 
ferromagnet [74,75]. Recall that the angular momentum (the 'charge') conjugate 
to the phase angle is the spin angular momentum of the overturned spins that 
form the skyrmion. In the quantum disordered 'insulating' phase, each skyrmion 
has a well defined integer-valued 'charge' (number of overturned spins) much 
like we found when we quantized the U{1) zero mode for the plane angle of 
a single isolated skyrmion in eq. (1.276). There is an excitation gap separating 
the energies of the discrete quantized values of the spin. 

The 'superfluid' state with broken U{1) symmetry is a totally new kind of 
spin state unique to non-colinear magnets [74, 75] . Here the phase angle is 

^^In fact in the quantum path integral this term has no effect except for time histories in 
which a 'vortex' encircles site j causing the phase to wind (pj{h/3) = ipjiO) ±27r. We explicitly 
ignore this possibility when we maJje the harmonic approximation. 
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wcll-dcfincd and the number of overturned spins is uncertain. The off-diagonal 
long-range order of a superfluid becomes 

{b]h) {e'^^'e-"^") (1.291) 

or in the spin language^* 

(s+{f)s-{f')). (1.292) 

Thus in a sense we can interpret a spin flip interaction between an electron 

and a nucleus as creating a boson in the superfluid. But this boson has a 
finite probability of 'disappearing' into the superfluid 'condensate' and hence 
the system does not have to pay the Zeeman price to create the flipped spin. 
That is, the superfluid state has an uncertain number of flipped spins (even 
though S'tQ^. commutes with H) and so the Zeeman energy cost is uncertain. 

In classical language the skyrmions locally have finite (slowly varying) x 
and y spin components which act as effective magnetic fields aroimd which the 
nuclear spins precess and which thus cause to change with time. The key 
here is that s'^ and can, because of the broken U{1) symmetry, fluctuate very 
slowly (i.e. at MHz frequencies that the nuclei can follow rather than just the 
very high Zeeman precession frequency). 

Detailed numerical calculations [75] show that the Skyrme lattice is very 
efficient at relaxing the nuclei and 1/Ti and is enhanced by a factor of ^ 10^ 
over the corresponding rate at zero magnetic fleld. We expect this qualitative 
distinction to survive even above the Skryme lattice melting temperature for 
the reasons discussed earlier. 

Because the nuclear relaxation rate increases by orders of magnitude, the 
equilibration time at low temperatures drops from hours to seconds. This means 
that the nuclei come into thermal equilibrium with the electrons and hence the 
lattice. The nuclei therefore have a well-defined temperature and contribute to 
the specific heat. Because the temperature is much greater than the nuclear 
Zeeman energy scale A ^ 1 mK, each nucleus contributes only a tiny amount 
~ /cg^ to the specific heat. On the other hand, the electronic specific heat 
per particle ~ fce -r^ is low and the electron density is low. In fact there are 

^ fermi 

about 10® nuclei per quantum well electron and the nuclei actually enhance the 
specific heat more than 5 orders of magnitude [67]! 

Surprisingly, at around 30 mK there is a further enhancement of the specific 
heat by an additional order of magnitude. This may be a signal of the Skyrme 
lattice melting transition [67,75,80], although the situation is somewhat murky 
at the present time. The peak can not possibly be due to the tiny amount of 
entropy change in the Skyrme lattice itself. Rather it is due to the nuclei in the 
thick AlAs barrier between the quantum wells. 

^^There is a slight complication here. Because the XY spin configuration of the skyrmion 
has a vortex-like structure (s+) = (s^ -fis") winds in phase around the skyrmion so the 'bose 
condensation' is not at zero wave vector. 

^^For somewhat complicated reasons it may be that the barrier nuclei are efficiently dipole 
coupled to the nuclei in the quantum wells (and therefore in thermal equilibrium) only due to 
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1.12 Double-Layer Quantum Hall Ferromagnets 
1.12.1 Introduction 

We learned in our study of quantum Hall ferromagnets that the Coulomb in- 
teraction plays an important role at Landau level filling factor v = 1 because 
it causes the electron spins to spontaneously align ferromagnetically and this 
in turn profoundly alters the charge excitation spectrum by producing a gap.^° 
A closely related eff'ect occurs in double-layer systems in which layer index is 
analogous to spin [43,44,81]. Building on our knowledge of the dynamics of 
ferromagnets developed in the last section, we will use this analogy to explore 
the rich physics of double-layer systems. 

Novel fractional quantum Hall effects due to correlations [82] in multicom- 
ponent systems were anticipated in early work by Halperin [42] and the now 
extensive literature has been reviewed in [43]. There have also been recent in- 
teresting studies of systems in which the spin and layer degrees of freedom are 
coupled in novel ways [83,84]. 

As described in this volume by Shayegan [45] , modern MBE techniques make 
it possible to produce double-layer (and multi-layer) two-dimensional electron 
gas systems of extremely low disorder and high mobility. As illustrated schemat- 
ically in Fig. (1.34), these systems consist of a pair of 2D electron gases separated 
by a distance d so small {d ^ lOOA) as to be comparable to the typical spacing 
between electrons in the same layer. A second type of system has also recently 
been developed to a high degree of perfection [85]. These systems consist of 
single wide quantum wells in which strong mixing of the two lowest electric sub- 
bands allows the electrons to localize themselves on opposites sides of the well 
to reduce their correlation energy. We will take the point of view that these 
systems can also be approximately viewed as double-well systems with some 
effective layer separation and tunnel barrier height. 

As we have already learned, correlations are especially important in the 
strong magnetic field regime because all electrons can be accommodated within 
the lowest Landau level and execute cyclotron orbits with a common kinetic 
energy. The fractional quantum Hall effect occurs when the system has a gap 
for making charged excitations, i.e. when the system is incompressible. Theory 
has predicted [42,82,86] that at some Landau level filling factors, gaps occur in 
double-layer systems only if interlayer interactions are sufficiently strong. These 
theoretical predictions have been confirmed [87]. More recently work from sev- 
eral different points of view [88 93] has suggested that inter-layer correlations 
can also lead to unusual broken symmetry states with a novel kind of sponta- 
neous phase coherence between layers which are isolated from each other except 
for inter-layer Coulomb interactions. It is this spontaneous interlayer phase co- 
herence which is responsible [43,51,73,94] for a variety of novel features seen in 

the critical slowing down of the electronic motion in the vicinity of the Skyrme lattice melting 
transition. 

^"Because the charged excitations are skyrmions, this gap is not as large as naive estimates 
would suggest, but it is still finite as long as the spin stiffness is finite. 
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Figure 1.34: Schematic conduction band edge profile for a double-layer two- 
dimensional electron gas system. Typical widths and separations are W d 
100 A and are comparable to the spacing between electrons within each inversion 
layer. 

the experimental data to be discussed below [44,81]. 
1.12.2 Pseudospin Analogy 

We will make the simplifying assumption that the Zccman energy is large enough 
that fluctuations of the (true) spin order can be ignored, leaving out the pos- 
sibility of mixed spin and pseudospin correlations [83,84]. We will limit our 
attention to the lowest electric subband of each quantum well (or cquivalently, 
the two lowest bands of a single wide well). Hence we have a two-state system 
that can be labeled by a pseudospin 1/2 degree of freedom. Pseudospin up 
means that the electron is in the (lowest electric subband of the) upper layer 
and pseudospin down means that the electron is in the (lowest electric subband 
of the) lower layer. 

Just as in our study of ferromagnetism we will consider states with total 
filling factor v = v-^ + = 1. A state exhibiting interlayer phase coherence 
and having the pseudospins ferromagnetically aligned in the direction defined 
by polar angle 9 and azimuthal angle (p can be written in the Landau gauge just 

as for ordinary spin 

W = n {cos(^/2)4^ + sin(0/2)e^^4} |0). (1.293) 

k 

Every k state contains one electron and hence this state has v = 1 as desired. 
Note however that the layer index for each electron is uncertain. The amplitude 
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to find a particular electron in the upper layer is cos(6'/2) and the amplitude to 
find it in the lower layer is sin(6'/2)e*'''. Even if the two layers are completely 
independent with no tunneling between them, quantum mechanics allows for the 
somewhat peculiar possibility that we are uncertain which layer the electron is 
in. 

For the case of ordinary spin we found that the Coulomb interaction pro- 
duced an exchange energy which strongly favored having the spins locally paral- 
lel. Using the fact that the Coulomb interaction is completely spin independent 
(it is only the Pauli principle that indirectly induces the ferromagnetism) we 
wrote down the spin rotation invariant effective theory in eq. (1.224). Here we 
do not have full SU(2) invariance because the interaction between electrons in 
the same layer is clearly stronger than the interaction between electrons in op- 
posite layers. Thus for example, if all the electrons are in the upper (or lower) 
layer, the system will look like a charged capacitor and have higher energy than 
if the layer occupancies are equal. Hence to leading order in gradients we expect 
the effective action to be modified slightly 



The spin stiffness Ps represents the SU(2) invariant part of the exchange energy 
and is therefore somewhat smaller than the value computed in eq. (1.231). The 
coefficient /? is a measure of the capacitive charging energy. The analog of the 
Zeeman energy A represents an external electric field applied along the MBE 
growth direction which unbalances the charge densities in the two layers. The 
coefficient t represents the amplitude for the electrons to tunnel between the 
two layers. It prefers the pseudospin to be aligned in the x direction because 
this corresponds to the spinor 



which represents the symmetric (i.e. bonding) linear combination of the two well 
states. The state with the pseudospin pointing in the —x direction represents the 
antisymmetric (i.e. antibonding) linear combination which is higher in energy. 

For the moment we will assume that both t and A vanish, leaving only the 
(3 term which breaks the pseudospin rotational symmetry. The case (3 < 
would represent 'Ising anisotropy'. Clearly the physically realistic case for the 
capacitive energy gives /3 > which represents so-called 'easy plane anisotropy.' 
The energy is minimized when = so that the order parameter lies in the 

•'^Wc have taken the charging energy to be a local quantity characterized by a fixed, wave 
vector independent capacitance. This is appropriate only if m''{f* ) represents the local charge 
imbalance between the layers coarse-grained over a scale larger than the layer separation. 
Any wave vector dependence of the capacitance will be represented by higher derivative terms 
which we will ignore. 



£= - 




(1.294) 




(1.295) 
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XY plane giving equal charge densities in the two layers. Thus we arc left with 
an effective XY model which should exhibit long-range off-diagonal order'^^ 



The order is 'off-diagonal' because it corresponds microscopically to an operator 



which is not diagonal in the basis, much as in a superfluid where the field 
operator changes the particle number and yet it condenses and acquires a finite 
expectation value. 

One other comment worth making at this point is that eq. (1.297) shows 

that, unlike the order parameter in a superconductor or superfluid, this one 
corresponds to a charge neutral operator. Hence it will be able to condense 
despite the strong magnetic field (which fills charged condensates with vortices 
and generally destroys the order). 

In the next subsection we review the experimental evidence that long-range 
XY correlations exist and that as a result, the system exhibits excitations which 
are highly collective in nature. After that we will return to further analysis and 
interpretation of the effective Lagrangian in eq. (1.294) to understand those 
excitations. 

1.12.3 Experimental Background 

As illustrated by the dashed lines in fig. (1.34), the lowest energy eigenstates 
split into symmetric and antisymmetric combinations separated by an energy 
gap AsAS = 2t which can, depending on the sample, vary from essentially zero 
to many hundreds of Kelvins. The splitting can therefore be much less than 
or greater than the interlayer interaction energy scale, = e^/ed. Thus it 
is possible to make systems which are in either the weak or strong correlation 
limits. 

When the layers are widely separated, there will be no correlations between 
them and we expect no dissipationless quantum Hall state since each layer has 
[95] v = 1/2. For smaller separations, it is observed experimentally that there 
is an excitation gap and a quantized Hall plateau [81,85,96]. This has either 
a trivial or a highly non-trivial explanation, depending on the ratio Asas/-£'c- 
For large Asas the electrons tunnel back and forth so rapidly that it is as if 
there is only a single quantum well. The tunnel splitting Asas is then analogous 
to the electric subband splitting in a (wide) single well. All symmetric states 
are occupied and all antisymmetric states are empty and we simply have the 
ordinary i/ = I integer Hall effect. Correlations are irrelevant in this limit and 
the excitation gap is close to the single-particle gap Asas (or hcuc, whichever 
is smaller). What is highly non-trivial about this system is the fact that the 

^■^At finite temperatures ^'(f') will vanish but will have long-range algebraically decaying 
correlations. Above the Kosterlitz-Thouless phase transition temperature, the correlations 
will fall ofi' exponentially. 



^f(f) = {rrf{r) +im'^{f)). 



(1.296) 



^-(f) = (s+(r)) = (^|(f)^4(f)) 



(1.297) 
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Figure 1.35: Phase diagram for the double layer QHE system (after Murphy et 
al. [81]). Only samples whose parameters lie below the dashed line exhibit a 
quantized Hall plateau and excitation gap. 



V = \ quantum Hall plateau survives even when AgAs ^ E^. In this limit 
the excitation gap has clearly changed to become highly collective in nature 
since the observed [81,85] gap can be on the scale of 20K even when Asas ~ 
1 K. Because of the spontaneously broken XY symmetry [51,73,88,89,92], 
the excitation gap actually survives the limit Asas — ^ 0! This cross-over from 
single-particle to collective gap is quite analogous to that for spin polarized single 
layers. There the excitation gap survives the limit of zero Zeeman splitting so 
long as the Coulomb interaction makes the spin stiffness non-zero. This effect 
in double-layer systems is visible in fig. (1.35) which shows the QHE phase 
diagram obtained by Murphy et al. [44,81] as a function of layer-separation and 
tunneling energy. A v ~ I quantum Hall plateau and gap is observed in the 
regime below the dashed line. Notice that far to the right, the single particle 
tunneling energy dominates over the coulomb energy and we have essentially 
a one-body integer QHE state. However the QHE survives all the way into 
Asas = provided that the layer separation is below a critical value d/fe ~ 2. 
In this limit there is no tunneling and the gap is purely many-body in origin and, 
as we will show, is associated with the remarkable 'pseudospin ferromagnetic' 
quantum state exhibiting spontaneous interlayer phase coherence. 

A second indication of the highly collective nature of the excitations can be 
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Figure 1.36: The charge activation energy gap, A, as a function of tilt angle 

in a weakly tunneling double- layer sample (Asas = 0.8K). The solid circles 
are for filling v = 1, open triangles for v = 2/3. The arrow indicates the 
critical angle Oc- The solid line is a guide to the eye. The dashed line refers 
to a simple estimate of the renormalization of the tunneling amplitude by the 
parallel magnetic field. Relative to the actual decrease, this one-body effect is 
very weak and we have neglected it. Inset: Arrhenius plot of dissipation. The 
low temperature activation energy is A = 8.66K and yet the gap collapses at 
a much lower temperature scale of about 0.4K (1/T f« 2.5). (After Murphy et 
al. [81]). 




seen in the Arrhenius plots of thermally activated dissipation [81] shown in the 
inset of fig. (1.36) The low temperature activation energy A is, as already noted, 
much larger than Asas- If A were nevertheless somehow a single-particle gap, 
one would expect the Arrhenius law to be valid up to temperatures of order 
A. Instead one observes a fairly sharp leveling off in the dissipation as the 
temperature increases past values as low as ^ 0.05A. This is consistent with 
the notion of a thermally induced collapse of the order that had been producing 
the collective gap. 

The third significant feature of the experimental data pointing to a highly- 
ordered collective state is the strong response of the system to relatively weak 
magnetic fields B|| applied in the plane of the 2D electron gases. In fig. (1.36) 
we see that the charge activation gap drops dramatically as the magnetic field 
is tilted (keeping B± constant). 

Within a model that neglects higher electric subbands, we can treat the 
electron gases as strictly two-dimensional. This is important since can affect 
the system only if there are processes that carry electrons around closed loops 
containing flux. A prototypical such process is illustrated in fig. (1.37). An 
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Figure 1.37: A process in a double-layer two-dimensional electron gas system 
which encloses flux from the parallel component of the magnetic fleld. One 
interpretation of this process is that an electron tunnels from the upper layer 
to the lower layer (near the left end of the figure). The resulting particle- hole 
pair then travels coherently to the right and is annihilated by a subsequent 
tunneling event in the reverse direction. The quantum amplitude for such paths 
is sensitive to the parallel component of the field. 
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electron tunnels from one layer to the other at point A, and travels to point B. 
Then it (or another indistinguishable electron) tunnels back and returns to the 
starting point. The parallel field contributes to the quantum amplitude for this 
process (in the 2D gas limit) a gauge-invariant Aharonov-Bohm phase factor 
exp (27ri$/$o) where <& is the enclosed flux and $o is the quantum of flux. 

Such loop paths evidently contribute significantly to correlations in the sys- 
tem since the activation energy gap is observed to decrease very rapidly with 
falhng by factors of order two or more until a critical field, iJy ~ 0.8T, 
is reached at which the gap essentially ceases changing [81]. To understand 
how remarkably small B*^ is, consider the following. We can define a length L[j 
from the size of the loop needed to enclose one quantum of flux: _L||_Bj| d = $o- 
(L||[A] = 4.137 X 10yd[k]B^^[T].) For B* = 0.8T and d = 150A, = 2700A 
which is approximately twenty times the spacing between electrons in a given 
layer and thirty times larger than the quantized cyclotron orbit radius £ = 
{hc/eB^y^"^ within an individual layer. Significant drops in the excitation gap 
are already seen at fields of O.IT implying enormous phase coherent correlation 
lengths must exist. Again this shows the highly-collective long-range nature of 
the ordering in this system. 

In the next subsection we shall briefly outline a detailed model which explains 
all these observed effects. 

1.12.4 Interlayer Phase Coherence 

The essential physics of spontaneous inter-layer phase coherence can be exam- 
ined from a microscopic point of view [51,73,90 -92] or a macroscopic Chern- 
Simons field theory point of view [51,73,88,89], but it is perhaps most easily 
visualized in the simple variational wave function which places the spins purely 
in the XY plane [51] 



Note for example, that ii ip = then we have precisely the non-interacting 
single Slater determinant ground state in which electrons are in the symmetric 
state which, as discussed previously in the analysis of the effective Lagrangian 
in eq. (1.294), minimizes the tunneling energy. This means that the system has 
a definite total number of particles [v = 1 exactly) but an indefinite number 
of particles in each layer. In the absence of inter-layer tunneling, the particle 
number in each layer is a good quantum number. Hence this wave function 
represents a state of spontaneously broken symmetry [51,88,89] in the same 
sense that the BCS state for a superconductor has indefinite (total) particle 
number but a definite phase relationship between states of different particle 
number. 

In the absence of tunneling {t = 0) the energy can not depend on the phase 
angle tp and the system exhibits a global U{1) symmetry associated with con- 
servation of particle number in each layer [88]. One can imagine allowing tp to 




(1.298) 
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vary slowly with position to produce excited states. Given the U{1) symmetry, 
the effective Hartree-Fock energy functional for these states is restricted to have 
the leading form 



The origin of the finite 'spin stiffness' ps is the loss of exchange energy which 
occurs when (p varies with position. Imagine that two particles approach each 
other. They are in a linear superposition of states in each of the layers (even 
though there is no tunneling!). If they are characterized by the same phase 
if, then the wave function is symmetric under pseudospin exchange and so the 
spatial wave function is antisymmetric and must vanish as the particles approach 
each other. This lowers the Coulomb energy. If a phase gradient exists then 
there is a larger amplitude for the particles to be near each other and hence the 
energy is higher. This loss of exchange energy is the source of the finite spin 
stiffness and is what causes the system to spontaneously 'magnetize'. 

We see immediately that the U{1) symmetry leads to eq. (1.299) which de- 
fines an effective XY model which will contain vortex excitations which interact 
logarithmically. [97, 98] In a superconducting film the vortices interact loga- 
rithmically because of the kinetic energy cost of the supercurrents circulating 
around the vortex centers. Here the same logarithm appears, but it is due to 
the potential energy cost (loss of exchange) associated with the phase gradients 
(circulating pseudo-spin currents). 

Hartrcc-Fock estimates [51] indicate; that the spin stiffness ps and hence 
the Kosterlitz-Thouless (KT) critical temperature are on the scale of 0.5 K in 
typical samples. Vortices in the (p field are reminiscent of Laughlin's fractionally 
charged quasiparticles but in this case carry charges zt^e and can be left- or 
right-handed for a total of four 'fiavors' [51,73]. It is also possible to show [51,94] 
that the presence of spontaneous magnetization due to the finite spin stiffness 
means that the charge excitation gap is finite (even though the tunnel splitting 
is zero). Thus the QHE survives [51] the limit AgAS — ^ 0. 

Since the 'charge' conjugate to the phase </? is the z component of the pseudo 
spin S^, the pseudospin 'supercurrent' 



represents oppositely directed charge currents in each layer. Below the KT tran- 
sition temperature, such current flow will be dissipationless (in linear response) 
just as in an ordinary superfluid. Likewise there will be a linearly dispersing col- 
lective Goldstone mode as in a superfiuid [51,73,88-90] rather than a mode with 
quadratic dispersion as in the SU(2) symmetric ferromagnet. (This is somewhat 
akin to the difference between an ideal bose gas and a repulsively interacting 
bose gas.) 

If found, this Kosterlitz-Thouless transition would be the first example of 
a finite-temperature phase transition in a QHE system. The transition itself 
has not yet been observed due to the tunneling amplitude t being significant in 
samples having the layers close enough together to have strong correlations. As 




(1.299) 



(1.300) 
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we have seen above however, significant effects which imply the existence of long- 
range XY order correlations have been found. Whether or not an appropriate 
sample can be constructed to observe the phase transition is an open question 
at this point. 

Exercise 1.23 Following the method used to derive eg. (1.230), show that 
the collective mode for the Lagrangian in eq. (1.294) has linear rather than 
quadratic dispersion due to the presence of the /3 term. (Assume A = t = 
0.) Hint: Consider small fluctuations of the magnetization away from m = 
(1, 0, 0) and choose an appropriate gauge for A for this circumstance. 
Present a qualitative argument that layer imbalance caused by A does not fun- 
damentally change any of the results described in this section but rather simply 
renormalizes quantities like the collective mode velocity. That is, explain why 
the u = 1 QHE state is robust against charge imbalance. (This is an impor- 
tant signature of the underlying physics. Certain other interlayer correlated 
states (such as the one at total filling v = 1/2) are quite sensitive to charge 
imbalance [43].) 



1.12.5 Interlayer Tunneling and Tilted Field Effects 

As mentioned earlier, a finite tunneling amplitude t between the layers breaks 
the U{1) symmetry 



iJeff = / d^r 



i/9s|V(/?p — ntcosif 



(1.301) 



by giving a preference to symmetric tunneling states. This can be seen from the 
tunneling Hamiltonian 

HT = -tj d\ {v|(r)Vi(r) + VIWVtW} (1-302) 

which can be written in the spin representation as 

Ht = -2t I d^rS^'ir). (1.303) 



(Recall that the eigenstates of are symmetric and antisymmetric combina- 
tions of up and down.) 

As the separation d increases, a critical point d* is reached at which the 
magnetization vanishes and the ordered phase is destroyed by quantum fluctu- 
ations [51,73]. This is illustrated in fig. (1.35). For finite tunneling t, we will 
see below that the collective mode becomes massive and quantum fluctuations 
will be less severe. Hence the phase boundary in fig. (1.35) curves upward with 
increasing Asas- 

The introduction of finite tunneling amplitude destroys the U(l) symmetry 
and makes the simple vortex-pair configuration extremely expensive. To lower 
the energy the system distorts the spin deviations into a domain wall or 'string' 
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Figure 1.38: Meron pair connected by a domain wall. Each meron carries a 
charge e/2 which tries to repel the other one. 
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connecting the vortex cores as shown in fig. (1.38). The spins arc oriented in 
the X direction everywhere except in the central domain waU region where they 
tumble rapidly through 27r. The domain wall has a fixed energy per unit length 
and so the vortices arc now confined by a linear 'string tension' rather than 
logarithmically. We can estimate the string tension by examining the energy 
of a domain wall of infinite length. The optimal form for a domain wall lying 
along the y axis is given by 



ip{'r) = 2 arcsin [tanh {Xx)], 
where the characteristic width of the string is 



t 



The resulting string tension is 

tPs 



To = 8 



27r^2 



(1.304) 



(1.305) 



(1.306) 



Provided the string is long enough {RX 1), the total energy of a segment of 
length R will be well-approximated by the expression 



E' 



pair 



4i? 



ToR. 



(1.307) 



This is minimized at R* = -^c'^/ATq. The linear confinement brings the charged 
vortices closer together and rapidly increases the Coulomb energy. In the limit 
of very large tunneling, the meron pair shrinks and the single-particle excitation 
(hole or extra spin-reversed electron) limit must be recovered. 

The presence of parallel field B|| field can be conveniently described with the 
gauge choice 



An = xBuZ 



(1.308) 



where z is the growth direction. In this gauge the tunneling amplitude trans- 
forms to 



t t e'*^" 
and the energy becomes 
1 



H 



I 



cos [if — Qx) 



(1.309) 



(1.310) 



where Q = 2tt/L\\ and Ly is the length associated with one quantum of flux 
for the loops shown in fig. 1.37. This is the so-called Pokrovsky-Talopov model 
which exhibits a commensurate-incommensurate phase transition. At low By, 
Q is small and the low energy state has ip w Qx; i.e. the local spin orientation 
'tumbles'. In contrast, at large Bu the gradient cost is too large and we have 
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(fi ss constant. It is possible to show [51,94] that this phase transition semiquan- 
titatively explains the rapid drop and subsequent leveling off of the activation 
energy vs. seen in fig. (1.36). 

Exercise 1.24 Derive eq. (1.304) for the form of the 'soliton' that minimizes 
the energy cost for the Hamiltonian in eq. (1.301). 
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Appendix A 

Lowest Landau Level 
Projection 



A convenient formulation of quantum mechanics within the subspace of the 
lowest Landau level (LLL) was developed by Girvin and Jach [26], and was 
exploited by Girvin, MacDonald and Platzman in the magneto-roton theory of 
collective excitations of the incompressible states responsible for the fractional 
quantum Hall effect [29] . Here we briefly review this formalism. See also Rcf. [8] . 

We first consider the one-body case and choose the symmetric gauge. The 
single-particle eigenfunctions of kinetic energy and angular momentum in the 
LLL are given in Eq. (1.76) 

where m is a non-negative integer, and z = {x + iy)/i. From (A.l) it is clear 
that any wave function in the LLL can be written in the form 

,Piz) = f{z)e-^ (A.2) 

where f{z) is an analytic function of z, so the subspace in the LLL is isomorphic 
to the Hilbert space of analytic functions [8,26,99]. Following Bargman [26,99], 
we define the inner product of two analytic functions as 

if,g) = ldt,{z)r{z)g{z), (A.3) 

where 

diJ.{z) = (27r)"^ dxdy e — -. (A.4) 

Now we can define bosonic ladder operators that connect (j)m to (l>m±i (and 
which act on the polynomial part of (j)m only): 



a 



z 

71' 



(A.5a) 
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V2 



dz' 



(A.5b) 



so that 



9) 

if, a 9) 



\/m + 1 ipm+i, 
(a^ f,9)- 



(A.6a) 

(A.6b) 
(A.6c) 
(A.6d) 



All operators that have non-zero matrix elements only within the LLL can be 
expressed in terms of a and . It is essential to notice that the adjoint of 
is not z* /\/2 but a = ^/2d/dz, because z* connects states in the LLL to higher 
Landau levels. Actually a is the projection of z*/^/2 onto the LLL as seen 
clearly in the following expression: 

(/,i^5) = (^ /,5) = (a^ f,9) = {f,ag). 
So we find 

^=2^, (A.7) 

where the overbar indicates projection onto the LLL. Since z* and z do not 
commute, we need to be very careful to properly order the operators before 
projection. A little thought shows that in order to project an operator which is 
a combination of z* and z, we must first normal order all the z*^s to the left of 
the z's, and then replace z* by z*. With this rule in mind and (A.7), we can 
easily project onto the LLL any operator that involves space coordinates only. 
For example, the one-body density operator in momentum space is 

n — ^ p-»qT — p-iiq'z+qz') — ^ p-^i^* p-ii'^ 
where A is the area of the system, and q = Qx + iQy Hence 

where 

Tg = e-'«^-5«*^ (A.9) 
is a unitary operator satisfying the closed Lie algebra 

r^Tfe = r,+fce^«^^ (A.lOa) 

\Tq,Tk\ = 2irg+feSin— — , (A.lOb) 
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z X (-q) 




(^ = 2% 




A 



z X q 



Figure A.l: Illustration of magnetic translations and phase factors. When an 
electron travels around a parallelogram (generated by TqTkT-qT-k) it picks up a 
phase (j) = 2t:-^ = q A k, where $ is the flux enclosed in the parallelogram and 
$0 is the flux quantum. 



where g A A; = £^(q x k) ■ z. We also have TqTk T-qT-k = e*'^*^. This is a 
familiar feature of the group of translations in a magnetic field, because qAk is 
exactly the phase generated by the flux in the parallelogram generated by q£'^ 
and \s.£^. Hence the r's form a representation of the magnetic translation group 
[see Fig. (A.l)]. In fact Tq translates the particle a distance x q. This means 
that different wave vector components of the charge density do not commute. 
It is from here that non-trivial dynamics arises even though the kinetic energy 
is totally quenched in the LLL subspace. 

This formalism is readily generalized to the case of many particles with spin, 
as we will show next. In a system with area A and N particles the projected 
charge and spin density operators are 



where Tq{i) is the magnetic translation operator for the ith particle and S^ is the 
lith component of the spin operator for the ith particle. We immediately find 
that unlike the unprojected operators, the projected spin and charge density 




(A.lla) 




(A.llb) 
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operators do not commute: 

[pk,S!;] = e J sin j ^ 0. (A.12) 

This implies that within the LLL, the dynamics of spin and charge are entangled, 
i.e., when you rotate spin, charge gets moved. As a consequence of that, spin 
textures carry charge as discussed in the text. 



Appendix B 

Berry's Phase and 
Adiabatic Transport 



Consider a quantum system with a Hamiltonian which depends on a set of 
externally c;ontrolled parameters represented by the vector R. Assume that for 
some domain of R there is always a finite excitation gap separating the ground 
state energy from the rest of the spectrum of H^. Consider now the situation 
where the parameters R{t) are slowly varied around a closed loop in parameter 
space in a time interval T 

R{0) = R{T). (B.l) 

If the circuit is transversed sufficiently slowly so that h/T ^ Amin where Amin 
is the minimum excitation gap along the circuit, then the state will evolve 
adiabatically. That is, the state will always be the local ground state \E' of the 

H(t) 

instantaneous Hamiltonian H^^^-^ . Given the complete set of energy eigenstates 
for a given R 

Hs'i'^i'' ^ e'^i^¥i\ (B.2) 

R R R R ^ ^ ' 

the solution of the time-dependent Schrodinger equation 

iJ^^^=Hn^t)HrM) (B.3) 

is 

rC\t) 

The adiabatic approximation consists of neglecting the admixture of excited 
states represented by the second term. In the limit of extremely slow variation 
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of R{t), this becomes exact as long as the excitation gap remains finite. The 
only unknown at this point is the Berry Phase [49] ^{t) which can be found 
by requiring that ■^(r, t) satisfy the time-dependent Schrodinger equation. The 
LHS of eq. (B.3) is 



dtp{r,t) 
'^—Qi- 



^ (r) 



if we neglect the aj{t) for j > 0. The RHS of eq. (B.3) is 
within the same approximation. Now using the completeness relation 



e-;., tpif,t) 



(B.5) 



(B.6) 



i=o 



^R^' R 



(B.7) 



In the adiabatic limit we can neglect the excited state contributions so eq. (B.5) 
becomes 



(0) 



dRf" m 



R{t) 



V'- 



This matches eq. (B.6) provided 



(B.8 



(B.9) 



The constraint 



^(0) 



1 guarantees that 7 is purely real. 



Notice that there is a kind of gauge freedom here. For each R we have a 
different set of basis states and we are free to choose their phases independently. 
We can think of this as a gauge choice in the parameter space. Hence 7 and 7 
are 'gauge dependent' quantities. It is often possible to choose a gauge in which 
7 vanishes. The key insight of Berry [49] however was that this is not always 
the case. For sonic; problems involving a closed-circuit F in parameter space the 
gauge invariant phase 



TBerry 



dt 7 



dR" 



d 



dRt' 



(B.IO) 



is non-zero. This is a gauge invariant quantity because the system returns to 
its starting point in parameter space and the arbitrary phase choice drops out 
of the answer. This is precisely analogous to the result in electrodynamics that 
the line integral of the vector potential around a closed loop is gauge invariant. 
In fact it is useful to define the 'Berry connection' A on the parameter space by 



A''{R) = i 



d 
dRM 



(B.ll) 
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which gives the suggestive formula 

^Berry = f dR -Air). (B.12) 

Notice that the Berry's phase is a purely geometric object independent of the 
particular velocity R^(t) and dependent solely on the path taken in parameter 
space. It is often easiest to evaluate this expression using Stokes theorem since 
the curl of ^ is a gauge invariant quantity. 

As a simple example [49] let us consider the Aharonov-Bohm effect where 
A will turn out to literally be the electromagnetic vector potential. Let there 
be an infinitely long solenoid running along the z axis. Consider a particle with 
charge q trapped inside a box by a potential V 

The position of the box is moved along a closed path R{t) which encircles 
the solenoid but keeps the particle outside the region of magnetic flux. Let 

X^"' — R{t)^ he the adiabatic wave fimction in the absence of the vector 

potential. Because the particle only sees the vector potential in a region where 
it has no curl, the exact wave function in the presence of A is readily constructed 

where the precise choice of integration path is immaterial since it is interior 
to the box where A has no curl. It is straightforward to verify that '^'^^] 

R(t) 

exactly solves the Schrodinger equation for the Hamiltonian in eq. (B.13) in the 
adiabatic limit. 

The arbitrary decision to start the line integral in eq. (B.14) at R constitutes 
a gauge choice in parameter space for the Berry connection. Using eq. (B.ll) 
the Berry connection is easily found to be 

A''{R) = +^ Af'iR) (B.15) 
nc 

and the Berry phase for the circuit around the flux tube is simply the Aharonov- 
Bohm phase 

7Berry = ^ dR^" A^ = 27r^ (B.16) 

where $ is the flux in the solenoid and $0 = hc/q is the flux quantum. 

As a second example [49] let us consider a quantum spin with Hamiltonian 

H = -A{t) ■ S. (B.17) 

The gap to the first excited state is 7i|A| and so the circuit in parameter space 
must avoid the origin A = where the spectrum has a degeneracy. Clearly the 
adiabatic ground state has 



(0) 



<Sf^°A = hS-^. (B.18) 

^ / lAI ^ ^ 
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If the orientation of A is defined by polar angle and azimuthal angle (p, the 
same must be true for (S). An appropriate set of states obeying this for the 
case S = h is 



since these obey 

(^e,^ \S'\tl>0,^) = nS (^cos'^ I - sin^ = nS cose (B.20) 

and 

{iPe,v\S'' + iSy\^e,v) = (V-e.y \S+\ tpe,^) = hS sinO e"^. (B.21) 

Consider the Berry's phase for the case where A rotates slowly about the z axis 
at constant 6 

d_ 



iBerry ^ j 



/■27r 
' / 

Jo 



dip I cos - sm - e 
dip (1 — cos^) 



= -'I 

Jo 

/■27r /■! 

= -S d>p dcose' = -Sn (B.22) 
^0 ^cos e 

where ft is the solid angle subtended by the path as viewed from the origin of 
the parameter space. This is precisely the Aharonov-Bohm phase one expects 
for a charge ~S particle traveling on the surface of a unit sphere surrounding 
a magnetic monopole. It turns out that it is the degeneracy in the spectrum at 
the origin which produces the monopole [49] . 

Notice that there is a singularity in the connection at the 'south pole' = tt. 
This can be viewed as the Dirac string (solenoid containing one quantum of 
flux) that is attached to the monopole. If we had chosen the basis 

e-"' \^sj (B.23) 

the singularity would have been at the north pole. The reader is directed to 
Berry's original paper [49] for further details. 

In order to correctly reproduce the Berry phase in a path integral for the 
spin whose Hamiltonian is given by eq. (B.17), the Lagrangian must be 

C = hS {-m^y4^ + A^m^ + A(m^m^ - 1)} (B.24) 

where m is the spin coordinate on a unit sphere, A enforces the length constraint, 
and 

VmXA = m (B.25) 

is the monopole vector potential. As discussed in the text in section 1.10, this 
Lagrangian correctly reproduces the spin precession equations of motion. 
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